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ABSTRACT 

We present 24 /im and 70 /im MIPS (Multiband Imaging Photometer for 
Spitzer) observations of 70 A through M-type dwarfs with estimated ages from 
8 Myr to 1.1 Gyr, as part of a Spitzer guaranteed time program, including a 
re-analysis of some previously published source photometry. Our sample is se- 
lected from stars with common youth indicators such as lithium abundance, X- 
ray activity, chromospheric activity, and rapid rotation. We compare our MIPS 
observations to empirically derived Ks-[24] colors as a function of the stellar ef- 
fective temperature to identify 24 /im and 70/im excesses. We place constraints 
or upper limits on dust temperatures and fractional infrared luminosities with a 
simple blackbody dust model. 

We confirm the previously published 70 /im excesses for HD 92945, HD 
112429, and AU Mic, and provide updated flux density measurements for these 
sources. We present the discovery of 70 /im excesses for five stars: HD 7590, 
HD 10008, HD 59967, HD 73350, and HD 135599. HD 135599 is also a known 
Spitzer IRS (Infra-Red Spectrograph) excess source, and we confirm the excess 
at 24 /im. We also present the detection of 24 /im excesses for ten stars: HD 
10008, GJ 3400A, HD 73350, HD 112429, HD 123998, HD 175742, AT Mic, BO 
Mic, HD 358623 and Gl 907.1. We find that large 70 /im excesses are less com- 
mon around stars with effective temperatures of less than 5000 K (3.7^]''^%) than 
around stars with effective temperatures between 5000 K and 6000 K (21.4^5;7%), 
despite the cooler stars having a younger median age in our sample (12 Myr vs. 
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340 Myr). We find that the previously reported excess for TWA 13 A at 70 fim is 
due to a nearby background galaxy, and the previously reported excess for HD 
177724 is due to saturation of the near-infrared photometry used to predict the 
mid-infrared stellar fiux contribution. 

In the Appendix, we present an updated analysis of dust grain removal time- 
scales due to grain-grain colhsions and radiation pressure, Poynting-Robertson 
drag, stellar wind drag and planet-dust dynamical interaction. We find that drag 
forces can be important for disk dynamics relative to grain-grain collisions for 
L/ii'/L*<10~^, and that stellar wind drag is more important than P-R drag for 
K and M dwarfs, and possibly for young (<1 Gyr) G dwarfs as well. 

Subject headings: circumstellar matter planetary systems: formation 



1. INTRODUCTION 



Nearly twenty-five years ago, the Infrared Astronomical Satellite (IRAS) launched the 
study of in frared excesses around stars that we attribute to mature extra-solar planetary sys- 



tems fe.g..lRhee et al.ll2007:IZuckerman 



2001 



Faiardo-Acosta et al.ll2000l:lMannings fc Barlow 



1998: Backman fc Paresc 



1993 



Walker fc Wolstencroft 



1988 



Aumann et al.lll984j ). Parent 



bodies in a planetary system - ~1 m and larger aggregates of rock and ice analogous to 
asteroids and Kuiper-belt objects in our own Solar System - collide and produce dusty 
debris that is heated by incident stellar radiation. This optically thin "debris disk" re- 
emits the absorbed radiation at infrared wavelengths; it is detecte d in excess of the ex - 
pected stellar radiation around ~15% of main sequen ce stars (e.g., Lagra.nge et aPboOO ). 



Ground-based infrared and sub-millimet er efforts (e.g.. iLestrade et al 



2006 



Liu et al 



2005 



Laureiis et al. 



20021: 



Habing et al. 



2004 



Weiriberger et al.| 20041: Song et al. l2002h, the Infrared Space Observ atory (ISO; de Muizon 



200 ll : ISpangler et al.l l200ll ) and the Spitzer Space 



Telescope ( IWerner et al.ll2004j ) have discovered >100 nearby stars with infrared excess. These 





RebuU et al. 


2008 


2005a 


). Large samples 



RebuU et al.l l2008l : iMever et al.l l2007l : ISmith et al.l l2006l : ILow et all l2005l : IChen et al. 



useful in determining overall trends such as disk frequency, disk infrared luminosity, and 



Carpenter et al.ll2009l: '. 


iillenbrand et al.ll2008 


:lTrilling et a 


.II2OO8: 


Wvatt 20081: ICautier et al 




20081; 


WvattI 2007 


{ Beichman et al. 2006al,bl: 


Brvden et al. 


200d: 


Su et al. 


2006|; 


Rieke et al 




2005 


Chen et al.ll20058 


) and in clusters and associations fe.g., Currie et al. 


2OO9I; 


Dahm & Carpenter 


2009 


Cieza et al. 


200^ 


: Currie et al. 2008bl: iHernandez et al. 2008: Cieza & Bahbei 2007 
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Currie et al. 


2007; 


Dahm & Hillenbrand 


2007; 


Gorlova et al. 


2007; 


Hernandez et al. 


2007a 


b; 


Carpenter et al.l 


2006; 


Chen et al. 


l2005b 


Low et al. 


2005; 


Stauffer et al. 


2005; 


Gorlova et al. 


2004 


). High-contrast, high-resolution direct and coronagraphic imaging from a number of 



telescopes has spatially resolved ~15 debris disks around solar-type stars in scattered and/or 
thermal e mission that show disk structures and clearings indicative of po ssible unseen plan 



Kalas et al. 


2004 


Stapelfeldt et al. 


2004; 


Holland et al.ll2003 


Wahhai et al.l 


2000; 


Schneider et al.l 


19991; IJavawardhana et al. 


19981; 


Koerner et al. 


1998^ 


ported discoveries of Fomalhaut b ( 


Kalas et al. 


I2OO8 


Ichiane et al. 


2OO9I 



2005 



Heap et al. 



12003 



The recent re- 
and a possible 



planet orbiting Beta Pic (jLagrange et al.ll2009l ) support the hypothesis that planets can di- 
rectly influence debris disk structure. Spatially resolved systems provide further information 
about the dust dynamics and evolution, but discerning overall trends is limited by the small 
number of resolved systems. 

In parallel with infrared observations and high-contrast imaging of young stars, over 
300 extrasolar planets have been discovered, primarily with the radial velocity technique, in 
the past decade, but also through direct imaging, transits, microlensing, and pulsar timing 



(iMarois et al. 



200 



Beaulieu et al. 



Kalas et al. 



20081; iButler et al.lbood : iMarcv et al l bood ; bond et al.l[2004 



2006 ; Wolszczan &: Frail Il992l . and references therein). The synergy of these 



planet-finding methods with the study of debris disks can yield new insights into the overall 
circumstellar architecture of exo-planetary systems. For example, radial velocity discovered 
extrasolar planets have revealed trends such as the Jovian planet frequency - host star 
metallicity correlation (iGonzalez 1119971 ). This correlation is not observed for debris disks 
( iBeichman et al.ll2006al ). 



Fundamental questions remain about the dust dynamics, properties, and evolution of 
debris disks around young stars. For example, the evolutionary time-scales for both primor- 
dial and debris disks appear to be dependent on spectral type. Optically thick primordial 
disks are common (~50%) at ages of ~1 M yr around stars with spectral types of A through 



M ( iHaisch et al.l 1200 ll ; iMeyer et al.l Il997l ) . These primordial disks provide the material to 
form planets and parent bo dies that in turn can ge nerate secondary debris disks. In the 3-5 
Myr Upper Sco association, ICarpenter et al.l (|2006l ) observes that primordial disks continue 
to persist in an optically thick state around K and M dwarfs, but have already transi- 
tioned to optically thin disks around earlier type stars. It might follow that debris disks 
should similarly persist for longer around K and M dwarfs due to lower stellar luminosi- 
ties and masses ( ^A.2). Debris disks can persist for stars older than ~1 Gyr around A,F , 



and G-type stars fiTrilling et al.ll2008l ; iMever et al.ll2007l : ISu et al.ll2006l : iBrvden et al.ll2006 
Rieke et al.ll2005l ). However, for K and M stars older than ~50 Myr, there i s a relative 



paucity of stars known with infrared excess associated with debris disks (e.g., iRhee et al. 
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20071 : iPlavchan et al.ll2005l ). Before Spitzer, selection effects from a flux-limited survey could 



partially account for this apparent difference in disk frequency and disk lifetime; debris disks 
around K and M dwarfs are thought to be rel atively colder and fainter than arou nd solar- 
type stars, due to the lower stellar luminosities (IRhee et al.ll2007l : IZuckerman Il200ll ). Indeed, 
two c older debris disks have been reported around M dwarfs with sub-millimeter observa- 
tions (ILestrade et al.ll2006l : iLiu et al.l l2004j ). More sensitive observations with the Spitzer 
Space Telescope and ground-based efforts continue to identify a relati ve lack of bright 70 



nm excesses (F7n/F^^,2-10) around K and M stars older than ~50 Myr (ITrilling et al.ll2008 



Gautier et al.ll2008l : iBeichman et al.ll2006al ). This paper addresses and confirms this dearth 



of relatively bright excesses at 70 fim around stars cooler than our Su n, but find tha t smal l 
(~10-15%) flux excesses at 24 fim are common at ages of < 100 Myr. iForbrich et al.l (120081 ) 
also identifies 9 new candidate M dwarf debris disks with large 24 /xm excesses in the ~40 
Myr open cluster NGC 2547, perhaps probing t he epoch of icy planet formation for M dwarfs 
JCurrie et al.lboOSbl : iKenvon fc Bromlevll2008l l 

Large ensemble infrared studies of young stars remain a valuable probe to address re- 
maining questions about the dynamics, properties, and evolution of debris disks. In this 
paper, we present the results of a Spitzer guaranteed time program to search for debris 
disks. We have obtained MIPS 24 fim and 70 /im observations of 70 young, nearby A-M 
stars with ages rangin g from 8 Myr t o ^1 G yr. Preliminary res u lts for a subset of these stars 
have been reported in lRebuU et al.l (120081 ): I Chen et al.l (j2005al ): ILow et al.l (120051 ). In §2 we 
present our sample and selection criteria. In §3 we present the observations and our data 
analysis procedures. In §4, we present our analysis to evaluate the significance of our detec- 
tions and non-detections, and place constraints on the disk properties of observed stars with 
and without infrared excess. In §5 we present our results on stars with identified infrared 
excesses. In §6, we discuss individual stars with excess in our sample, and the transition 
from grain-grain collisions to drag-force dominated disks. In §7, we present our conclusions. 
In the Appendix we present an updated analysis of the different time-scales for dust grain 
removal in debris disks from radiative forces, stellar wind forces, and grain-grain collisions. 
We also consider the importance of planet-dust interaction relative to these processes. 



SAMPLE SELECTION 



Our sample of 70 stars were selected as being relatively young based upon a variety of 
youth indicators. For F-M stars, selection criteria include high X-ray activity from ROSAT 
detection s, high lithium abundance. Calciu m II H&K line chromospheric activity, and rapid 
rotation (lGaidoslll998l : iGaidos et al.l 120001 ). For A stars, we select candidates based on in- 
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ferred youth from Stroemgren photometry and dis tances inferred from Hipparcos parallaxes 
( jSong et al.ll200ll : ISongll200d : IStaufferllpriv.comm.1 ). with the exception of HD 123998 which 



was selected for its apparent IRAS exc ess (confirmed in this work t o be a false-positive, po- 
tentially due to IRAS beam-confusion; ISong||2000l : ISong et al.ll200ll ). We select 27 stars with 
stellar effective temperatures less than 5000 K, 28 stars with stellar effective temperatures 
between 5000 K and 6000 K, and 15 stars with stellar effective temperatures greater than 
6000 K. We use these sub-groups within the sample throughout our analysis. The choice of 
5000 K and 6000 K to delineate the sub-groups is motivated by the different sample selection 
techniques applied in each temperature regime, the X-ray activity, and estimated ages (§4.4). 
The sub-samples also highlight differences in our results between solar- type (5000-6000 K), 
low-mass (< 5000 K) and high- mass radiative atmosphere stars (> 6000 K). Many of the 
stars in our sample are members or candidate members of young moving groups and associ- 
ations. In Table 1 we list our targets along with literature spectral types, distances, cluster 
and moving group memberships, and other youth indicators. 



OBSERVATIONS 



We obtained 24 /im and 70 fim o bservations of our sample using the Multiband Imag- 



i ng Photometer for Spitzer (MIPS, iRieke et al.l |200J) on the Spitzer Space Telescope 



( iWerner et al.ll2004l ) in photometry mode (default scale). Each of our targets was observed 
in 2004 or 2005, using integration times of 48.2 s at 24 ^, m. and 186-465 s at 70 /xm. The data 
were reduced using the MIPS Data Analysis Tool (DAT. lGordon et al.ll2005l ). which performs 
the necessary steps need to produce a final mosaic of the target. We also applied additional 
reduction techniques to r emove instrurnental artifacts and flux calibration as described by 
Engelbracht et al.l (120071 ): iGordon et al.l (120071 ) at 24 and 70 um, respectively. We adopt flux 
density calibration factors at 24 /im and 70 yum of 0.0454 and 702 respectively to convert 
from instrumental units to MJy per steradian. We also reduce and combine the "bed" image 
data from the Spitzer Science Center "S13"/M0PEX/APEX pipeline. The "S13" pipeline 
processes and removes image artifacts such as bad-pixels, cosmic rays, latent images, verti- 
cal array column "jailbar s" from calibrating "stim flashes" for the 70 /im detector, and the 
"droop effect" at 24 /im (ICordon et aDl2005h . 



We extract photometry using Point-Spread Function (PSF) fits to the co-added image 
data with the DAT. We examine the 70 /zm images for any cirrus contamination. We compare 



PSF- d erived flux densities to aperture photometry derived with the DAT as in ISu et al. 



( 120061 ). as well as to flux densities derived from the MOPEX/APEX pipeline using both 
aperture photometry and fitting to empirically derived Pixel-Response Functions (PRFs, 
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Makovoz fc Marleaul l2005l ). In all four methods, we set a minimum detection threshold 
signal-to-noise ratio (S/N) of 3. For both sets of comparison aperture photometry, we use 
a circ ular aper t ure w ith radii of 14. "94 (6 pixels) at 24 fim and 29. "5 (3 pixels) at 70 /xm 



as m iSu et al.l (120061 ) . and use scalar aperture corrections of 1.143 and 1.694 respec tively 
inferred from Spitzer Tiny Tim models of the point-spread function (IKrist et al.l 120021 ). We 
estimate background noise from the standard deviation of flux density in background annuli 
with inner and outer radii of 29". 55 and 42".33 (12-17 pixels) and 39".40 and 68". 95 (4-7 
pixels) at 24 yum and 70 fim respectively. This includes detector noise and noise due to cirrus 
structures present in the image. If the observations were taken in regions with high cirrus, 
then the detection limits are dominated by cirrus noise. Our sample cirrus levels span from 
very low (< 5 MJy/sr at 70 fim) to high (e.g., ~20 MJy/sr at 70 fim for HD 125158). 

In Table 2, we tabulate adopted 24 /xm and 70 fim photometry, propagated uncertainties, 
S/N, expected photospheric flux densities (§4), and references in which Spitzer photometry 
for these sources has been previously published. All photometry in Table 2 assumes a A^^ 
flux density across the bandpass. We find good agreement in the derived photometry for 
all four techniques at 24 /xm, and for S/N ^ 5 at 70 fim. At 70 /xm, neither the DAT 
nor MOPEX/APEX demonstrate greater sensitivity at low S/N. As noted in Table 2, there 
are seven instances where one or more methods do report a detection at 70 /im at S/N>3 
whereas the other methods do not. These results at 70 fim are below the adopted threshold 
for the detection of excess (§4.2) in all seven cases, and do not affect our results. For 
the 24 /im excess source AT Mic, a partially resolved ~3" binary, we use a DAT aperture 
photometry rather than PSF-fitting at 24 fim to measure the combined flux density from 
both components. 

At both 24 and 70 /im we use the DAT PSF photometry when available due to the 
smaller flux density calibration uncertainties. At 24 /im, we detect all of our sample targets. 
At 70 /im, we detect 19 of our 70 sample targets with S/N>3 at the expected positions in- 
ferred from the 24 /im source position, unless otherwise noted in Table 2. For non-detections 
at 70 /im, we calculate 3-a upper limits by measuring the standard deviation of flux den- 
sity within a circular aperture with a radius of 29.5" at the exp e cted source locatio r i. Ou r 
photometry are consistent with preliminary results in I Chen et al.l (l2005al ): iLow et al.l (120051 ) . 
but we obtain better overall precision and sensitivity. 



ANALYSIS 



In j^ 4.1, we use a minira ization procedure to flt PHOENIX NextGen model photo- 



spheres ( iHauschildt et al.lll999al ) to optical and near-infrared photometry to estimate stellar 
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effective temperatu re and radii. In j^4.2 , we compare the Ks-[24] colors to the relation em- 
pirically derived by iGautier et al.l (120081 ) from nearby (<5pc) disk-less main-sequence G, K 
and M dwarfs. In §4.3, we place constraints on the circumstellar material for detected and 
non-detected mid-infrared excesses. In §4.4 we estimate ages, masses, and calculate X-ray 
luminosities for our sample. 



4.1. Modeling of Stellar Photospheres 



We obtain photo-electric optical UBVRcIc photometry from the literature, with B and 
V magnitudes primarily from the Tyc ho catalog, tran s formed to t he Johnson photometric 
system using the corrected relations in lMamajek et al.l (120061 . |2002| ). We ob tain JHK., near 



infra red photometry from the 2MASS All-Sky Catalog of Point Sources (ISkrutskie et al. 



20061 ) . In several cases, our sources saturate the 2MASS detectors, and fo r two of these 
sources we obtain J and K band magnitudes from lMorel fc MagnenatI (119781 ). Due to satu- 
ration, we do not include 2MASS J-band photometry for HD 112429, HD 16555, HD 125158, 
and HD 144197, and H-band photometry for HD 112429. We also do not include optical 
photometry for several sources when the literature magnitudes appear inconsistent with the 
magnitudes in other bands - Rc-band magnitudes for HD 10008, TWA 18, and HK Aqr, and 
I-band magnitudes for SAO 145139, TWA 12, and HK Aqr. We tabulate and cite all the 
photometry in Table 3. 

PHOENIX NextGen models are used to fit the observed optical and near-IR photometry 
rather than a blackbody, since the stellar SEDs are very dif fe rent from that of a blackbody at 
optic al and near-IR wavelengths ( Hauschildt et al. 1999a| jbl: lOray et aL 1992 : MuUan et al. 



19891 ). We assu me solar metallicity as is typical for main-seq uence stars in the solar neigh- 
borhood (e.g., ISoone et allbood : brilling fc Landoltl boooh . We assume log g=4.5 as is 
typical for main sequence stars, and Av=0. Our sample includes some K and M dwarfs 
that have not yet reached the main sequence and will have lower surface gravities, but this 
has a negligible effect on our derived stellar radii and effective temperatures. To fit the 
NextGen spectra to the available photometry, we compute a minimization as a function 
of effective temperature and normalization. We integrate model photospheres as a function 
of wavelength across effective bandpasses to compare to observed photometry for a given 
band. We assign a weight of 1 to each band, and we normalize to K^. 

For effective temperatures (T*) greater than 3600 K, we interpolate (linearly in T^) 
the NextGen model photospheres (provided in increments of 200 K) in increments of 20 K 
to improve the fits. We find excellent agreement with effective temperatures infer red from 
literature spectral types (lLuhmanlll999l : iLuhman fc Riekelll998l : iHartigan et al. I ll994h . For T* 
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less th an 3600 K, NextGen photospheres lack near-infrared opacity sources (iHauschildt et al. 



1999al ). As a result, model photosphere fits for these M dwarfs tend towards lower effective 



temperatures (up to 200 K) than those implied from their previously identified spectral 
types. The resulting fits in turn overpredict the observed 24 /xm photometry by as much 
as 20%. For these sources, we fix the effective temperature in increments of 100 K to the 
literature spectral type. This approach provided the mo st consistent results in predicted 24 
lim photospheric flux density with the method in §4.2 (iGautier et al.ll2008l ). We tabulate 
the derived stellar parameters from our model photosphere fits in Table 3. 



4.2. K,-[24] and [24]-[70] Colors 



To identify excesses at 24 and 70 /im, we use the effective temperatures derived in §4.1 
and the K,,-[ 24] color relation deri ved from nearby (<5pc) disk-less main sequence G, K and 
M dwarfs in iGautier et al.l (120081 ) to predict the expected photopheric flux densities at 24 
fim: 



Ks-[2A] 



-0.2915 + [1.267 X 10"^ + (2.389 x 10"^ + (2.2 x 10"^ T,)^] ^ , < 5025K 
, T, > 5025/^ 

(1) 

where [24] is the magnitude converted from the 24 /xm photometry with a zero-point flux of 
7.14 Jy, and T^, is the effective stellar temperature in Kelvin. At 70 /im we use a Rayleigh- 
Jeans extrapolation from the expected stellar flux at 24 fim. 

We calculate the signiflcance of the deviation of the measured photometry from expected 
photospheric flux densities, x = (-F;/ (observed) — Fj, (predicted) )/(T,y at 24 and 70 fim, where 
cTj, includes the detection S/N and flux calibration uncertainties. We list these values in 
Table 2 and plot them in Figure 1. We identify sources with X70>5 as 70 fim excesses, X24>5 
as 24 /im excesses. In Figure 2 we plot Fjo/F^, as a function of effective temperature for 
our sample. Figure 2 includes several strong detections of 70 fim disks reported here for the 
flrst time. In Figure 3 we plot the observed Ks-[24] colors for our sample as a function of 
effective stellar temperature. Figure 3 includes eleven detections of 24 /im disks reported 
here for the flrst tim e. While some of the new 24 /xm disks have been previously analyzed in 



Chen et al.l (l2005al ). the improved calibration uncertainties at 24 fim with the DAT pipeline 



have improved the sensitivity to small levels (e.g. ~10-15%) of flux density excess. 



- 9 - 



4.3. Dust Disk Constraints 

We list the derived dust disk parameters and 70 /xm color corrected photometry in Table 
4. To estimate the physical disk parameters - fractional infrared excess, dust temperature, 
dust mass, and orbital distance - for the stars with and without excesses in our sample, 
we assume that the debris disks are optically thin. After subtracting off the predicted 
photospheric flux density, we fit a single-temperature dust blackbody model to the 24 fim 
and 70 /xm excess. The derived fractional infrared excesses are consistent with our model 
assumption of an optically thin disk. We estimate dust orbital radii assuming that the 
dust is in thermal equilibrium. We color correct the 70 fim photometry based on the dust 
temperature, using the color corrections outlined in the Spitzer Observer's Manual. The 
color corrections are <1% at 24 fim, and ~10% at 70 /xm. For the purposes of estimating 
upper limits to the observed fractional infrared excesses, we fit to 3-0" upper limits for non- 
detections at 70 /im. In Figure 4 we plot the model photosphere and dust blackbody fits for 
stars with 70 fim detections. For the sources with 24 /im excesses and 70 fim non-detections, 
we estimate the fractional infrared excess at 24 fim with the assumption that L//j~z/F,^(24 
yum), after subtracting from F^, the expected stellar flux density contribution. We list this 
fractional infrared excess in Table 4 in addition to the upper limit derived from the 70 fim 
non-detection. 



We estimate dust masses as in I Chen et al.l (j2005al ). inferring the average grain size 



(a) using a coUisional cascade grain size distribution power law exponent of -3.5, with a 
minimum grain size set by the radiative blowout radius a.b ((a) =1.32 a;,), and grain density 
of Pg = 2.5 g/cm~^ (Equation A22 in §A.2.3). We assume /5 = 1/2 is the threshold criterion 
for radiative blowout, where P is the ratio of radiation pressure to gravity for the purposes 
of calculating a?, (Equation A18, §A.2.2). For the stars in our sample, only significant stellar 
winds of ~ IOOOM0 will affect the blowout radius (§A.1,§A.2.2). S uch a large mass-loss 



rate is not supported by theoretical considerations and observations (IHolzwarth fc Jardine 



20071 ). and hence only grain blowout by radiation pressure is relevant (§A.l). We derive 
dust disk masses between 2.8x10^^ and 5.2x10"^ lunar masses for the disks in our sample. 
We estimate minimum parent body masses required to produce this dust between 3xl0~^ 
and 1.1 lunar masses (Equation A20). The required minimum parent body masses are ~10^ 
times the derived dust disk mass. These results imply both a broad dynamic range in mass 
of debris disks, and a diverse variety of asteroid or Kuiper Belt extra-solar analogs. 
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4.4. X-Ray Luminosity, Age, Mass, and Stellar Wind Mass-Loss Rate 

Estimates 

We estimate properties of the host stars in our sample that are relevant to the investi- 
gation of debris disk evolution and spectral type dependence. In Table 4, we summarize the 
derived ages, age determination methods, masses, L^/L*, and stellar wind mass-loss rates. 



4.4- i- X-Ray Luminosity 



We derive X- ray apparent flux den sities from ROSAT PSPC detected count-rates using 
the conversion of ISchmitt et al.l (Il995l ): 



Apparent Flux = # cnts/s x (8.31 + 5.3 x HRl) x 10"^^ ergs/cm^/s 



(2) 



where # cnts/s is the ROSAT count rate, and HRl is the ROSAT hardness ratio #1. 
We require that optical source positions fall within the 1-a ROSAT position uncertainties 
to associate the ROSAT detections with our sample stars. We infer the fractional X-ray 
luminosity hx/'^* and the X-ray flux density Fx at the stellar photosphere using the derived 
stellar properties (radius, luminosity) from synthetic spectral fltting in §4.1. 



4-4-2. Stellar Ages 



We use flve main methods to estimate ages for the stars in our sample - young moving 
group membership (<100 Myr), X-ray activity, Ca II H & K activity, rotation periods, 
and fltting to theoretical isochrones. First, we assign ages based on cluster and association 
memberships to AU Mic, AT Mic, and HD 358623 (Beta Pic moving group; 12 Myr), to 
TWA association sources (8 Myr), and AB Dor moving group members AB Dor, GJ 3400A 
and LO Peg (75±25 Myr). For sources with measured X-ray activity, rotation periods 
or Ca II H &: K activity, we estim ate ages using activity-age correlation relations derived 
Mamajek &: Hillenbrand! (120081 ). We use equations A3, 12-14, and 3 respectively, with 



m 



appropriate limits (e.g., B-V<1.2, -5 < Log(Lx/L,,) < -4.3). When two or three of these 
activity indexes are measured, we adopt the mean and standard deviation as the flnal age 
estimate and uncertainty. When only a single indicator is available, we adopt an uncertainty 
of 20% for Call H&K and X-ray activity. Fo r rotation-based age estimates, we propagate 
the uncertainties in the coefficients of the fit inlMamaiek fc Hillenbrandl (120081 ) to Equations 
12-14, the relation first presented in iBarned (120071 ). We obtain age estimates for 46 of our 
70 target stars using these first four age determination methods. 



For eight of the remaining stars with T* < 6000 K and X-Ray acti vity implying ages 



<300 Myr, we fit the derived stellar radii and temperature in §4.1 to ISeiss et al.l (120001 ) 
isochrones. We adopt average ages and uncertainties from the range of isochrone grid points 
that are consistent to within 5% of the derived stellar radius and temperature. The 5% 
criterion is chosen based on the estimated radius uncertainties in the fits in §4.1, the S/N 
of Hipparcos measured parallax, and to recover reasonable age uncertainties (~50%). If a 
source is an unresolved binary, we adjust the derived radius appropriately before fitting to 
isochrones to avoid systematically under-estimating the age. Three of the eight stars have 
measured Lithium abundances that support the estimated ages. For stars that are implied 
to be young from measured activity indicators, the activity-based ages are consistent with 



Sone 


(2000 


); 


Sone; et al. 


(2001 


); 


Stauffer 


(priv.comm. 



We are unable to estimate ages for two target stars using the above five methods. For 
the possible Pleiades moving group member RE 1816+541 lacking a trigono metric parallax. 



we as si gn an age of 1 15 M yr based upon the Pleiades age estimates from iMessina et al. 



( 2003 ): IStauffer et al.l (1l998l ). Finally, we obtain an age estimate of 1.1 Gyr from iNordstrom 
( j2004j ) for HD 17240 and adopt an age uncertainty of 50%. 



In Figure 5, we plot the age distribution for our sample as a function of effective stellar 
temperature. The K and M stars in our sample are primarily selected from young moving 
groups whereas the A-G type stars are primarily selected from other indicators. As a result 
of these selection criteria, the stars in our sample cluster in age as a function of spectral type 
in Figure 5. The T=i,< 5000 K low- mass K and M dwarfs are younger than the 5000-6000 
K G dwarfs and the T=, > 6000 K stars, with median ages of 12.25 Myr, 340 Myr, and 200 
Myr respectively. This age bias limits comparisons across spectral types in our sample. 



4-4-3- Stellar Wind Mass-Loss Rates 



We wish to consider stellar wind drag from stellar wind protons as a disk dissipation 
mechanism (§A.l). In order to evaluate this mechanism, we estimate the stellar wind mass- 
loss rates for stars in our sample. If Fx<8xl0^ ergs/s/cm^ at the stella r surface (^ st ars 
approximately older than 700 Myr), we can use the relation of I Wood et al.l (120021 . |2005| ) to 
estimate stellar mass-loss rates, M„„: 
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where Mq = 2 x 10 ^'^Mq/ji. Seven stars with X-ray detections meet this criteria in 
our sample. We caution against estimating mass-loss rat es from X-Ray act ivity for stars 
younger than ~700 Myr or with Fx>8xl0^ ergs/s/cm^ Jwood et aliboosh . HD 135599 
has an abnormally low X-ray luminosity given its estimated age of 214±17 Myr, and also 
possesses a strong 70 /im excess. 



4-4-4- Stellar Mass 



We estimate stellar mass e s for our sample using the estimated stellar ages, effective 
temperatures, and ISeiss et al.l (120001 ) isochrones. We use the estimated stellar mass in the 
calculation s of du st disk properties and time-scales. For a give n age and c orresp onding 



Seiss et al.l (120001 ) isochrone, we linearly interpolate between the ISeiss et al.l ( 120001 ) stellar 



masses as a function of stellar effective temperature. 



5. RESULTS 



5.1. Stars with Excess 



At 70 /im, we report the discovery of 70 /xm excesses (x70>5) for HD 7590, HD 10008, 
HD 59967, HD 73350, and HD 135599. \ ye confirm the p r eviously detected MIPS excesses 
for HD 92945, HD 112429 and AU Mic JLow et al.l[2005l : Ichen et al.ll2005ah . We confirm 
that the Spitzer In fra-Red Spectrogra ph (IRS) detected excess for HD 135599 extends to 
longer wavelengths (ILawler et al.l 120091 ). The model photospheres and empirical colors from 
§4.1 and §4.2 produce consistent values of X7o at 70 /xm. We do not confirm the previously 
published 70 fim excess for TWA 13A, and instead attribu te the 70 /im emission to a nearby 
background galaxy (~10", Figure 6). I Chen et al.l (l2005al ) previously reported an excess for 
HD 177724 at 24 /im and 70 /im that wasn't easily described by a simple dust blackbody 
model for the excess. We instead attribute the previously reported excess to the large (0.28- 
0.36 mag) errors in the saturated J, H and K^-band 2MASS photometry, and find that the 
MIPS photometry is consistent with the predicted photospheric fiux density. For all stars 
with excesses at 70 /im, we PSF subtract the detected sources to look for extended residuals. 
A U Mic and HP 92945 s how extended structures, and detailed analysis will be presented 



m 



Bryden et al.l (lin prep. I ). The 70 /tm residuals for AU Mic and HD 92945 appear aligned 



with the position angles of the disks resolved in scattere d optical and near-infrared light 
jKalas et al.lbooi iMetchev et al.l[2005l : iKrist et al.ll2005al jbl). albeit at low S/N. 



At 24 /im, we report the discovery of eleven sources with 24 /im excesses (x24>5). We 
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identify flux density excesses for HD 10008 (16.6±2.0%), GJ 3400A (13.9±3.3%), HD 73350 
(9.5±1.9%), HD 112429 (11.1±1.9%), HD 123998 (10.5±2.7%), HD 135599 (10.7±2.8%), HD 
175742 (14.5±2.2%), AT Mic (15.1±4.2%), BO Mic (14.6±3.3%), HD 358623 (11.4±2.3%) 
and Gl 907.1 (9.1±2.3%). We attribute tliese small excesses to warmer circumstellar material 
than responsible for a 70 /im-only excess. The above uncertainties in the flux density excess 
are slightly larger than implied by the X24 values in Table 2 because we include the uncer- 
tainties in the band flux density as an estimate of the model photosphere flux density 
uncertainty at 24 fim. Five of the 24 /xm excess sources have effective temperatures between 
3200 and 4800 K. Thus this is one of the first reports of debris disks around cool, low mass 
stars in addition to iForbrich et al.l (120081 ). 



5.1.1. Alternative Explanation for 24 j^m Excess 

We investigate the possibility that the small excesses we detect at 24 fim are not due 
to re-processed radiation from dusty debris disks. HD 10008, HD 73350, HD 112429 and 
HD 135599 also possess 70 /im excesses, supporting the assertion that our data and analysis 
are sensitive to ~10% excesses at 24 fim. For the seven sources with 24 /im excesses and 
without 70 fim detections, the upper-limits at 70 fim are consistent with a debris disk SED. 
Furthermore, excesses at 24 fim are more likely at younger ages. The stars with 24 fim 
excesses we identify are generally much younger than our overall sample age distribution. 

However, an unknown low-mass spectroscopic companion - e.g., a M5 companion to 
a KO star - could produce an apparent ~10% flux excess at 24 fim. The occurrence of 
spectroscopic binaries in the field with the mass ratio required to reproduce t he 24 fim 



exces s is much less than the 10% (7/70) of stars for which we observe this excess ( iHogeveen 



19921 ). Nonetheless, we do not rule out low-mass companions for these sources. BO Mic and 



GJ 3400A possess published IRAQ photometry , but the [8]- [24] colors are inconclusive as 



to which scenario is favored (jChen et al.ll2005al ). Positionally coincident background AGN 



and serendipitous asteroids could also account for the observed 70 and 24 fim excess, but 
are statistically unlikely for any given star (<1%). >0.1 mag stellar variability at K^-band 
could also produce an apparent 24 fim excess for non-contemporaneous observations. 

We note that the distribution of Ks-[24] photospheric colors for low- mass stars is less 
constrained than it is for warmer stars. We compare the formal errors derived from the 
24 fim S/N and K^-band uncertainty, and compare it to the empirical Ks-[24] distribution 
after subtracting off the predicted photospheric color. We derive an empirical distribution 
uncertainty of 0.03 mag from a Gaussian fit, which is consistent with our formal propagated 
errors. 9 of the 11 24 fim excess sources lie >3-a in excess above the fit to the empirical 
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distribution. The 10th, HD 112429, also possesses a 70 /iin excess, and the 11th at 2.5-cr, is 
12.5 Myr 4240 K GJ 907.1. A 24 /iin excess for GJ 907.1 is entirely plausible, but the excess 
for Gl 907.1 is the least robust detection of a 24/im excess that we include in our results. 
For a sample size of 70 stars, we might expect ~1 false-positive detection of excess from 
statistical fluctuations. 



5.2. Disk Fraction 



tion is smaller than the disk fr action for A stars (33 = b5% 



For our sample, we find an overall disk fraction of 21. 4l^'^% (15/70) . This disk frac 



Su et aLlliood), sl i ghtly higher 



than for young G stars (14^7^% ICarpenter eTaDbood : ImUenbrand et al Iboosi: iMever et all 



20051). and higher than the disk fraction fo r older, field G stars (12-15%. iTrilling et al 



2008 



Bryden et al.l l2006l : iLagrange et al.l 120001). We compute the 1-cr uncertainties in the disk 
fraction using binomial statistics as in lBurgasser et al.l (120031 ). We list these values in Table 
5 for a variety of sub- samples. 



5.2.1. Spectral Type 



around stars with effective temperatures of less than 5000 K (3.7^^^%) 
with effective temperatures between 5000 K and 6000 K (21.4+|^% 



At 70 /im, we can detect debris disks that are generally >100% brighter than the 
expected photosphere (Figure 2). For the 15 stars in our sample with T^,>6000 K, we find 
a disk fraction of 13.3^4^6^%. T his is smaller t han the disk fraction derived from analogous 
studies of A-type stars (33±5%. ISu et al.ll2006l ). but the results are consistent to within 1-a 
given the small number statistics. In our sample 70 /xm excesses are ~6 times less common 

than around stars 
These two disk fractions 
differ at the ~2-cr level. We observe this discrepancy even though the cooler stars have a 
younger median age in our sample (12.25 Myr vs. 340 Myr). At younger ages, debris disks 
might be expected to be brighter by factors of a few. This paucity of late-type star debris 
disks at 70 /im has been found by, and is consistent with, several other more detailed stu dies 
JmUenbrand et al.lboosi : iGautier et aPboOsI : iTrilhng et allboosl : iBeichman et al.ll2006al . and 
references therein). The interpretation of our result is complicated and limited by several 
factors and selection effects in our sample that we consider in more detail in §A.3. 

At 24 /im, we can detect much smaller (~10-15%) excesses relative to the photosphere 
(Figure 3). When we include the detections of excesses at 24 /im, the disk fraction dependence 
on spectral type disappears. We find an overall disk excess fraction of 22.2~^^'q% for stars with 
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effective temperatures of less than 5000 K, and 25.0]'^g3% for stars with effective temperatures 
between 5000 K and 6000 K. This result suggests that debris disks are just as common 
around cool stars as they are around hotter stars, but the dust disks around hotter stars 
may generally be more massive and brighter relative to the stellar photosphere. This result 
is consistent w i th th e frequency of ~40 Myr M dwarfs with 24 /xm excesses identified in 



Forbrich et all (120081 ). 



5. 2. 2. Binarity 



Our sample includes 24 known binary systems (Table 1), and constitutes a sample that 
has been fairly well studied for AO, wide separation and spectroscopic companions. The 
binary completeness of our sample, however, is undoubtedly <100%. The incompleteness 
of the known binaries in our sample limit our interpretation, but we can make a tentative 
comparison. We do not identify a statistically significant difference between the disk fractions 
for binaries and single stars. We derive a s ingle star disk fracti on of 22.7^50% (10/44), and 
a binary disk fraction of 19. 2154% (5/26). Trilling et al. (j2007 ) carried out an investigation 
of debris disks around A and F-type binaries, and concluded that binaries have a higher 
disk fraction on average than single stars. We do not confirm this resu lt when investigating 
lower mass stars, with the aforementioned caveats. iTrilling et al.l (120071 ) identified that tight 
(separations <3 AU) and wide separation (>50 AU) binaries have a higher disk fraction than 
intermediate separation binaries (3-50 AU). Our sample of binaries is evenly split between 
intermediate separation (12) and wide separation binaries (14), and we find 3 and 2 debris 
disks in each category respectively. 



5.3. Age 

In Figure 7, we plot the fractional infrared excess as a function of stellar age in our 
sample. We over-plot a t7^ power law consistent with coUisional evolution of debris disks 



)A.2.1, 



e.g. 



WyattI I2OO8I . I2OO7I : iRieke et al.ll2005l ). The data are not consistent with this 



power-law evolution. However, the youngest stars in our sample (e.g. <100 Myr) are dom- 
inated by the coolest stars (T^, < 5000 K) and vice-versa. Consequently, this bias masks 
any true underlying trend in L/jj/L* vs. age. We caution against drawing any significant 



24 and 70 /xm has been i nvestigated in other more detailed studies 



2009 



s (e.g.. 


Carpenter et al. 


Carpenter et al. ( 


2OO9I) 



and iTrilling et al. rj2008h both note that the 70 /xm age evolution for the disk frequency and 
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disk brightne ss for solar-type stars does not de crease with age as strongly as is observed for 
A-type stars JSu et aliboOQ : iRieke et al.lboosh . 



Five of the seven 24 /im excess sources without large 70 fim excesses are younger than 
^50 Myr. This age range is consistent with the epoc h of terrestrial and icy planet formation 



[e.g. 



Kenyon fc Bromley! l2008l : ICurrie et al.l l2008bl ) , and is much younger than the median 
age of our sample (200 Myr). The average age of stars with T* < 5000 K is 21 Myr and 
60 Myr for stars with and without detected infrared excess respectively (at either 24 or 70 
fim). For comparison, the averages are 312 Myr and 395 Myr respectively for 5000 K < T,,, 
< 6000 K stars. The stars with excesses tend towards younger ages, and this is consistent 



with a coUisional-cascade production of dust from parent bodies ( IWyattll2008l . and references 
therein) . 



5.4. Stellar Rotation Period 



In Figure 8, we plot the fractional infrared excess as a function of stellar rotation period 
for our sample for all stars with T^,<6000 K. We use measured stellar rotation periods when 
available in Table 1. If not, we estimate the rotation period from measured v sini rotational 
velocities in Table 1, estimated radii in Table 3, and we assume an average rotation axis 
inclination with respect to our line of site of 30°. For HD 92945, we a ssume the star mat ches 
the inclination of the debris disk, estimated to be 25° from edge on (IKrist et al.ll2005bl ). 



We identify a weak correlation between rotational periods and observed fractional in- 
frared excesses for stars with excess. We derive a Spearman ra nk correlation coeffi cient 
Ps=0.35. Our sample spans an age range older than the sample in iRebuU et al.l (120081 ). who 
report a similar weak correlation between rotation and excess for stars in the 12 Myr Beta 
Pic moving group . Our o bservations support the premise presented in iRebuU et al.l (120081 ) 
and ICurrie et al.l (l2008al ) that massive primordial disks evolve into massive debris disks. 
Massive primordial disks regulate the initial angular momentum evolution of stars before 
the primordial disks dissipate, and the persisting excess-period correlation for debris disks 
is an echo of the primordial disk rotation c orrela ti on. A correlat i on be tween rotation and 
excess was also presented in IStauffer et al.l (120071 ) . IStauffer et al.l (120051 ) identified 3 debris 
disks in the ~100 Myr Pleiades, and did not identify any anti-correlation between excess 
and rotational velocity. 



-17- 



5.5. Activity Indicators 



In our sample, stars with 24 fim excess appear well-mixed with non-excess stars for X- 
ray saturated stars with log(Lx/L^,)>-4. Rotation periods between ~5.5-7 days are favored 
for the detection of 70 /im excess around stars with log(Lx/L^,)<-4. We do not identify any 
significant correlations of fractional infrared excess with either the lithium abundance or the 
activity index ^'hk- We do however note in §6.1.3 and Figure 9 that the bright 70 excess 
stars are outliers when computing the ratio of the rotation period to both log(Lx/L=K) and 
^HK- Given our small sample size and broad age and spectral type ranges, this trend is 
intriguing but inconclusive and possibly spurious. 

For the seven stars with Fx<8xl0^ ergs /s/cm^, we can use the X-ray lurainosity as a 



proxy for stellar wind mass loss rate (§A.2, iPlavchan et al.ll2005l : IWood et al.ll2005l ). We 



find that the smallest mass-loss rate in this sub-sample is seen for the only star with an 
excess - HD 135599. The other 6 stars possess excess upper limits considerably smaller 
than the excess detected for H P 135599. At a f i xed ag e, such an anti-correlation is expected 
from the steady-state model of [F 



the younger Sco-Cen association (I Chen et al 



avchan et a 



Jj2005|) and has been tentatively observed in 



2005bl ). Since we only have seven stars in this 



sub-sample and only one detected excess, we cannot draw any co nclusions except to no te 
that these observations do not contradict the steady-state model of iPlavchan et al.l (120051 ). 



6. DISCUSSION 

In §6.1, we present notes about the stars in our sample with the brightest excesses. In 
§6.2, we discuss the impact of various grain removal mechanisms on debris disk evolution 
and for the debris disks in our sample. These mechanisms are presented in detail in the 
Appendix. 



6.1. Individual Star Notes 



6.1.1. AU Mic 



Since the discovery of a resolved disk around AU Mic in iKalas et al.l (12004f). much 
work has been done to characterize and mode l AU Mic 's debris disk (IFitzgerald et al. 



2OO7I : Istrubbe fc Chianejbood : lOuillen et al1l2007l : iKrist et al.l[2005al jbl: iMetchev et al 



20051). 



Our updated 70 /im flux density measurement should not significantly alter interpreta- 
tions. We update the application of the stellar wind drag steady-state model to AU Mic 
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m iPlavchan et al.l (120051 ) . The discovery in IWood et al.l (120051 ) that X-ray activity does 
not correlate with stellar wind mas s-loss rates for Fx>8xl0^ ergs/s/cm^ invalidates the 
application in iPlavchan et al.l (120051 ) of a steady-state stellar wind drag model to AU Mic; 
AU Mic possesses an X-ray surface flux density of ~1.8xl0^ ergs/s/cm^. We estimate the 
grain-grain collision time-scale within the "birth ring" a t 43AU (the location of the dust- 
producing ring of parent bodies. IStrubbe fc Chiangll2006l ) to be ~10'^ times shorter than the 
stellar wind drag time-scale for M,.,,, = lOMfT^, a dust d i sk mass of 10~^Mm, and a disk width 
of 4.3 AU (Equation A2 3, §A.2, IPitzgerald et al.ll2007l : iMetchev et aPbood : lliu et all [2004 : 



Augereau &: BeustI l2006l ) . The evolution of the AU Mic debris disk is currently dominated 
by grain-grain collisions and the presence of any approximately Jovian mass planets (§A.2). 



6.1.2. HD 135599 

For this KO star (~1/3L0), the inferred mass-loss rate implies that the drag time-scale 
from stellar wind is ~10 times shorter than P-R drag (Equations A8-10, §A.1,§A.2),. The 
stellar wind drag time-scale is ~0.4 times the grain-grain collision time-scale (Equations 
A14,A23) for the dust radius (12 AU) and dust mass (6.6 x 10"^ Mmoon) estimated in Table 
4 (§4.3). This implies that stellar wind drag is dynamically important in the evolution of 
the dust in this disk, and makes HD 135599 a canonical example to study the transition 
between disk evolution dominated by grain-grain collisions and drag-force dominated disk 
evolution (§6.2 & §A2.2.3). The Spitzer IRS spectra of HD 135599 are consistent with the 
"power-law" excess pre dicted for a drag-for ce dominated debris disk (e.g. Fi,(dust)ocz/", §6.2 



and §A.2 in this paper. iLawler et al.l 120091 ). A more accurate dust mass measurement from 
sub-mm observations and high-contrast imaging to obtain radial surface brightness profiles 
is warranted to clarify this interpretation further. 

A power-law extrapolation of the IRS excess fo r HD 135599 predicts a flux at 70 fim that 



is smaller than what is observed (ILawler et al.ll2009l ). This implies that the entire disk for HD 



135599 is not described by a single dominant dust removal mechanism. Due to the different 
power-law dependencies on orbital distance for drag forces and grain-grain collisions, and 
because the collisional time-scale depends on the local dust disk density, the evolution of 
dust at different orbital radii can be dictated by different dust removal mechanisms (§A.2). 
For HD 135599, it is plausible that there is a ring of parent bodies producing dust through 
mutual collisions where the dust collision time-scales are short relative to the stellar wind 
drag time-scale, but not signiflcantly shorter. This would produce the "excess excess" seen 
at 70 fim. Then, a fraction of the dust produced in the ring spirals inward under the action 
of stellar wind drag into smaller orbital radii where the collisional time-scales become longer 
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than the drag force time-scales. This would then explain the observed IRS power-law excess. 



6.1.3. HD 7590, HD 73350, HD 92945, and HD 135599: An Activity-Excess 

Anti- Correlation ? 



We consider the large fractional infrared excesses (F7o/F^,>10) observed for HD 7590, 
HD 73350, HD 92945, and HD 135599 (the "bright disk stars"). For HD 7590, HD 73350 
and HD 135599, the stellar ages inferred from X-ray activity and tend towards older 
ages, whereas ages inferred from rotation periods tend towards younger ageil]. This trend 
is not present for stars with no excess, and follows directly by comparing the ratio of the 
three activity indexes. Within our sample, HD 7590, HD 73350, and HD 135599 are outliers 
when computing the ratio of the rotation period to both Log(R'j:^^) and Log(Lx/L^,) (Figure 
9). Given the small size of our sample, this result is inconclusive and may be spurious. A 
more rigorous treatment with a larger debris disk sample of solar-type stars with measured 
rotation periods and activity indexes is warranted. 

If HD 7590, HD 73350 and HD 135599 are younger than we estimate from X-ray and 
activity, then the activity for these stars are abnormally low. This presents the possibility 
that stars with large excesses are more likely to be less active at a given age. This anti- 
correlation between activity and excess is consistent with the possibility that stellar wind drag 
is the dominant grain removal mechanism (§A.l). At a fixed age, stellar wind drag predicts 
that less active sta rs will have brighter disks due to the longer wind drag time-scales (§A.2, 
Chen et al.l l2005bl ). However, at these ages of <700 Myr, X-ray activ ity is not a reliable 
proxy for stellar wind mass-loss rate as discovered in (jWood et al.ll2005l ). Alternatively, this 
anti-correlation is also consistent with massive debris disks evolvi ng from massive prim ordial 
disks, which in turn regulated the host star angular momentum (jCurrie et al.ll2008al ). 



A second hypothesis is that the "bright disk stars" are undergoing an epoch of en- 
hanced planetesimal collisions akin to the Late Heavy Bombardnient w i thin our own Solar 
System when it w as ~700 Myr old Jxieine et ahlbooi Ichambersi [20041 : IComes et aP 12005 



Strom et all 120051 ). HD 10008, HD 59967, and HD 123998 also fall within the estimated age 
range of 200-500 Myr. 



^For HD 7590, we estimate ages of 496, 427, and 463 Myr from X-ray, rotation period and indicators 
respectively; for HD 73350 we estimate ages of 613, 358, and 569 Myr from X-ray, rotation period and 
^'hk indicators respectively; for HD 92945, we estimate ages of 311 and 278 Myr from X-ray and 
activity respectively; for HD 135599, we estimate an age of 214 Myr from its rotation period and association 
membership, while we ignore the X-ray estimated age of 1.1 Gyr. 
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6.2. The Transition Between Grain-Grain Collisions and Drag Forces in Disk 
Dynamical Evolution — Observational Predictions 

Dust is produced in debris disks via a coUisional cascade of dust grains that result 
from the col l ision of parent bodies. T he collisions of parent bodies decay with time (e.g., 



WyattI I2OO8I : iKenyon fc Bromley! 120051 . and references therein) . The evolution of the ob- 
servable debris disk depends on the how the dust is removed from the system, primarily 
thought to occur through drag forces causing dust to spiral into the star (drag force dom- 
inated), radiative blowout as the end stage of a collisional cascade of dust grains (collision 
dominated), and/or planet-dust interactions (§A.2). Drag-force dominated disk evolution 
predicts a constant vertical optical depth as a function of orbital radius. This produces a 
power-law infrared spectrum, where F,^(dust)ocz/" over a wavelength regime consistent with 
the orbital distances of the in-spiraling dust in thermal equilibrium, and a varies from ~ 
-1- -3 depending on the wavelength dependence of the grain emissivity (§A.2.2). This is in 
contrast to a dust excess consistent with a single-temperature blackbody, such as produced 
by a single ring of dust (§A.2.1). Collision-dominated disks are necessarily relatively bright 
(log(L/ij/L*)~-3- -4). Resol ved collision-dominated debris disks are typically consistent 



with the ring geo me trv fe.g.. iKalas et al 



2008 



consistent with a single-temperature blackbody (e.g.. Morales et al.ll2009l : IChen et al.ll2006l ). 



Fitzgerald et all 120071: iTelesco et all 12005 



Krist et all r2005al jbl: IXalas et allbooi Istapelfeldt et allbooi). and the IRS excess spectr a 



To evaluate the importance of different grain removal mechanisms in our observed debris 
disks, we calculate the ratio of the drag time-scales to the grain-grain collision time-scales 
using equation A23 (§A.2.3). This ratio is approximately equal to or less than 1 at the 
minimum dust orbital radii and dust disk masses inferred for HD 10008, GJ 3400A, HD 59967, 
HD 123998, HD 175742, AT Mic, BO Mic, HD 358623 and Gl 907.1 (hereafter the "faint disk 
stars" ). This ratio suggests that dust could be removed under that action of drag forces before 
they are blown out from radiation pressure and collisional grinding. The fractional infrared 
excesses range from 6-80x10^® for the "faint disk stars". These "faint disk stars" possess 
flux excesses that are smaller t han have been typically investigated previously with Spitzer 



2008 



at 70 Mm (e.g. log(L rB/L* )>-4.ICarpenter et al.ll2009l : iHillenbrand et al.ll2008l : iTrilling et al. 



Bryden et al.l l2006l : iBeichman et al.l l2006al ) . Consequently, with this sample we are 
probing a dust dynamics regime previously unexplored. From the ratio of collision and 
P-R/wind drag time-scales, we infer that the transition between coUisional-cascade disk 
evolution to P-R drag disk evolution could occur at a fractional infrared luminosity as high 
as LjrI L^, ~ 8 X 10~^ (^A . 2.3). This result is in disagreement with the conclusions of IWyatt 
(120051 ) : iDominik fc Deciru (120031 ) which predict a transition at L/ij/L^, ~ 4 x 10^^, more than 



two orders of magnitude fainter. 
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The Spitzer IRS power-law spectrum for HD 135599 is a canonical example for the 
transition of the dominant dust dynamical mechanisms in a debris disk, with Lj^/L^, = 
1.0 X 10~^ (§6.1.2). A detailed IRS study of debris disks with a range of fracti onal infrared 



excess woul d further test this h y pothesis, and t his h as be en carried out in iLawler et al. 



J2009h and [Morales et all J2009h . Ilawler et all J2009h and [Morales et all J2009h find IRS 
excesses that are consistent with a power-law for excesses between Lir/L^ ~ 3 x 10~^ and 
LjFt/L^: ~ 3 X 10~^. This fractional infrared excess range is coincident with the fractional 
infrared excesses of the "faint disk stars." As noted for HD 135599 in §6.1.2, the transition 
from grain collision-dominated disk evolution to drag force-dominated disk evolution is not a 
sharp transition. Most disks in transition between collision time-scales and drag time-scales 
for dust removal will have different evolutionary dynamics at different orbital radii. Rings of 
parent bodies will produce dust locally that will initially be subjec t to grain-grain collisions . 



but drag forces can act to fill in the interior parts of the disk (e.g., IStrubbe &: Chiang||2006l ). 



For HD 7590, HD 73350, HD 92945, HD 112429, HD 135599, and AU Mic, stellar 
winds of ~3-600 times solar are necessary for drag force time-scales to overcome grain-grain 
collisions in dictating the disk evolution. While this is feasible for the ages of these stars, 
an important observational test is infrared spectra. Spitzer IRS observations of AU Mic and 
HD 92945 do not sho w the power law behavior expe cted if drag forces are the dominant grain 



removal mechanism (jChen et al. 



priv.comm. 



20061 ) . This implies the stellar wind mass- loss 



rates are less than ~50 Mq, and the evolution of these disks are dictated by the grain-grain 
collision time-scales. 

Finally, we note that planets of sufficient mass can truncate the inward spiral of dust 
grains under the action of drag forces before the dust reaches the sublimation radius (Equa- 
tion A30 and §A.2.4) . The presence of a power-law excess can be used to place an upper 
limit to the mass of any planet orbiting within the orbital radii of observed dust. 



7. CONCLUSIONS 

Infrared photometry provides basic constraints on fractional fiux excess, dust temper- 
ature, and dust orbital radii in debris disks associated with young planetary systems. We 
identify new 70 /im excesses for HD 7590, HD 10008, HD 59967, HD 73350, and HD 135599, 
and constrain their dust disk properties. We confirm the excesses at 70 fim previously iden- 
tified for HD 92945, HD 112429, and AU Mic. We present the discovery of 24 /im excesses 
for HD 10008, GJ 3400A, HD 73350, HD 112429, HD 123998, HD 175742, AT Mic, BO Mic, 
HD 358623 and Gl 907.1, and confirm the 24 /xm excess for HD 135599. 
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Within our 70 star sample selected for youth indicators, we find that large 70 //m 

excesses (L£;/L^>10~^) associated with debris disks are more common around ~200-500 
Myr solar-type stars with T^, > 5000 K than for ~8-50 Myr stars with < 5000 K. 
We find that small 24 yum excesses (L(;/L^,<10~^) associated with debris disks arc equally 
common among young stars of all spectral types. The tentative correlations we present with 
activity indicators are speculative given the small and non-uniform sample, but suggest a 
more detailed investigation is warranted. 

From arguments presented in the Appendix and in §6.2, we predict that drag-forces 
are relevant compared to grain-grain collisions for the dynamical evolution of disks around 
G-M dwarfs with Ljji/L^<10^'^. We find that stellar wind drag is more important than P-R 
drag for removing dust around K and M dwarfs with these fractional infrared excesses, and 
possibly for G dwarfs as well. 
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A. Appendix 



In §A.l, we present a brief overview of stellar wind drag. In §A.2 we summarize 
and update equations for the grain lifetimes due to grain-grain collisions, stellar wind 
forces, radiative forces, and the presence of planets. This section can serve as a refer- 
ence to quickly evaluate the relative importance of these process es in debris disk evolutio n 

in 
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In §A.3, we present five factors that affect the interpretation of the debris disk fraction as a 
function of spectral type in §5.2.1. 



A.l. Stellar Wind Drag Primer 



Plavchan et al.l (120051 ) and lMinato et al.l (120061 ) propose that stellar wind drag acting on 



dust grains in a debris disk could account for the relative paucity of late-type stars older than 
~50 Myr with detectable infrared excess. For our Sun, protons from the solar wind produce 
a radial o utward pressure, ~3xl0^ times smaller in magnitude than the analogous radiation 
pressure ( iGustafsonI Il994j ). Radiation pressure also gives rise to the Poynting- Robertson 
drag force that is reduced in magnitude from the radial component by a factor of forb/c, 
the ratio of the dust grain orbital velocity to the speed of light. Stellar wind pressure also 
has an analogous drag term, reduced in magnitude from the radial component by a factor 
of Vorb/vsw, where Vsw is the stellar wind velocity. This factor is the aberration angle of the 
protons seen by a dust grain as it orbits a star. A dust grain preferentially forward scatters 
stellar wind protons in the direction it is traveling, causing it to lose angular momentum 
and slowly spiral inwards towards the star. Since stellar wind velocities are on the order of a 
few hundreds of km/s, it follows that Vorb/vsw^^O^ Vorb/c. While stellar wind radial pressure 
is effectively negligible compared to radiation pressure for dust dynamics around solar-type 
stars, the drag terms are instead comparable in magnitude. 

For the Su n, the strength o f stellar wind drag is ~l/3^'^ the strength of Poynting- 
Robertson drag (jGustafsoru Il994j ) . From §A.2, we find that stellar wind drag will be more 
important than Poynting-Robertson drag for stars with mass-loss rates comparable to the 
Sun and stellar luminosities less than 1/3 L0 - e.g. K and M dwarfs - and for stars with 
enhanced mass-loss rates - e.g. young stars. As an example, for a ~l/3 Lq K dwarf, a mass- 
loss rate of 5 time solar would lead to a stellar wind drag force that is 5 times stronger than 
P-R drag. This in turn shortens dust lifetimes by a factor of 5, and lowers the observable 
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fractional infrared excess compa red to a disk w i th on ly P-R drag considered. This follows 
directly from Equation AlO and iPlavchan et al. J2005h . 



Quantitatively evaluating the relevance of stellar wind drag in low-mass star debris disk 
evolution is complicated by the uncertainty in the strength and evolution of stellar winds 
from low mass stars. Recent observations confirm that stellar winds from M dwarfs in M 
dwarf - white dwarf binaries are ^2-10 tira es the mass-lo ss rate from the S un at ages of 



several Gyr (ISchmidt et al.ll2007l : lDebesll2006f ). Additionally, IWoodetalJ (j2005|) observe that 
stellar wind mass-loss rates correlate with X-ray activity, implying that young stars also have 
higher mass-loss rates and hence shorter dust grain lifetimes under the a ction of stellar wind 
drag. However, for stars younger than ~500-700 Myr, the correlation in IWood et al.l (120051 ) 
breaks down, and the X-ray luminosity cannot be used as a direct proxy for measuring 
the stellar wind mass-loss rate. Young late-type stars exhibit enhanced X-ray activity, but 
this is not a direct indication of a correspondingly large radial stellar wind mass-loss rate. 
Measuring stellar wind mass-loss rates and their evolution from young (~5-500 Myr), low 
mass stars remains an open problem. 

J2005h to infer that M dwarfs older than ~700 



We can use the relation of IWood et al. 
Myr, such as the M dwarfs observed in 



Gautier et al.l (120081 ) . possess mass loss rates per 



units surface area that are comparable to solar-type stars. For a dust disk comparable in 
mass to the zodiacal dust cloud in our Solar System, the stellar wind drag time-scale will be 
shorter than both the P-R drag and grain-grain collision time-scales, and the disk evolution 
will be dominated by stella r wind drag. This can account for the l ack of older M dwarf s 
with observed debris disks (IGautier et al.l l2008l : iMinato et al.l l2006l : iPlavchan et al.l |2005| ) . 
We explore the role of stellar wind drag relative to other grain removal processes in more 
detail 5A.2. 



A. 2. Dust Grain Removal time-scales 



A. 2.1. Grain- Grain Gollisions 



T he average tim e for a grain to collide with anoth er grain is estimated with the expres- 



sion m 



Minato et al.l (120061 ): iNaiita fc WiUiamd (120051 ) as: 
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where a is the grain radius, r is the dust orbital distance, pg is the dust grain density, 
is the stellar mass, is the dust mass, Mq is the solar mass, is the Earth's mass, 
and {tio? /Sz) = l/Q^^jj, the ratio of the effective grain-grain collision cross-section to the 
geometric cross-section {Qqq\\ ~ 1). 

Equation Al relies on two important assumptions. First, the grains are assumed to be 
sufficiently inclined in their orbits to see the full vertical optical depth of the disk at the 
orbital radius of the grain (the grain traverses the full optical depth of the disk twice as it 
orbits the parent star). For optically thin debris disks, the average grain inclination should 
approximate the ratio of the scale height of the dust disk to the orbital distance, h/r. 

For the second assumption, the dust is assumed to be evenly distributed over a disk 
surface area ttDq^j,. This second assumption is only valid when the entire dust mass 
is interior to the orbital distance under consideration, ie r = Dqut- In this sense, t^q\\ in 
Equation Al represents the time-scale for all grains at the disk outer radius Dqut of size a 
and density pg in an entire disk of mass Md to undergo a collision. 



With approximately a dozen bri ght, spatially res o 



ved disks, the observed dust ge 



ometries are primarily narrow rings ( 


Kalas et al. 


2008: 


Fitzeeralc 


et al 


2007 


: Krist et al 




2005a b: Telesco et al. 2005: Marsh et al. 2005: S 


itapelfeldt et al. 


2004: 


Kalas et al.l 2004 




Holland et al. 20031: Wahhai et al.ll2003 


Heap et al 


2000: 


Schneider et al.lll999: 


Javawardhana et al. 


1998h. Here we update Equation Al and do awav with the second assumption to avoid mis- 



application of coUisional time-scales as a function of orbital radius, and to make more realistic 
estimates. 

We define the radial surface density profile Sr(r) such that: 



Dr 



Din 



Sr(r) X r dr 



(A2) 



where Djj^r and Dqut are the inner and outer disk radii respectively, and we have assumed 
azimuthal symmetry. Note we have incorporated the constant from the azimuthal integration 
Then the general form for Equation Al becomes: 



into T,r(r 



^coU - 2 X 10 yr (^^^j J J J 1^]^^ J [^-^ x (soAU)^ ) 

(A3) 

From Equation A2, if we set Dj^ = 0, and T,r{r) = Sq, we get that T,r{r) = Sq = 
2Mj:,/ D'^^j,^ and we recover Equation Al. 
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We consider a generic radial surface density profile of: 

S,(r) = Sor° (A4) 

with inner and and outer disk radii Du^ and Dqut, which can be set appropriately depending 
on a. From Equation A2, we find 



Md/HDout/Din) a = -2 



and equation A3 becomes: 



^coll - 2 X 10 yr 7 j (^Jo^ 



where: 

f (J^\ \( Dout \'^+^ I Din \°+^ 

ry _ ) \a+2) \ \rMAU I \bQ AU J 

\ 2\n{DouT/DiM) 



,a = -2 



(A6) 
(A7) 



The gamma factor in Equations A6 and A7, encapsulating the disk geometry, is impor- 
tant for estimating grain lifetimes. For example, consider a ring of dust at radius Dring and 
width AD ring=0. ID ring- The collision time-scale estimate in Equation Al will differ from the 
actual collision time-scale by a factor of 7 ~2ADri„g=0.2 for a=-2 or 0, and overestimate 
the collision frequency for any of the dust that is scattered or dragged out of the ring. As 
a second example, we consider a = —1.5 with Din=0 AU and Dout=200A\J . The collision 
time as a function of orbital distance will scale as r^, and the collision time-scales will be 
7 =8 times longer. 



A. 2. 2. Radiative and Stellar Wind Drag 



A drag force removes angular-momentum from an orbiting dust grain, causing it to 
spiral into the star. The time for a dust grain to spiral into a star under the action of P-R 
drag is estimated in iMinato et al.l (120061 ) as: 



TpR = 2x 10^ yr 



50AU 



a 



Pa 



Tia 



-1 



(A^ 



^g cm V \Cph^ 

where (na'^/Cph) = l/Qrad, where Qrad is the radiative coupling coefficient, the ratio of the 
effective to geometric cross-sections. Similarly, the time for a dust grain to spiral into a star 



under the action of stellar wind drag is estimated in 



Minatoetal 



fl2006f ) as: 



Tsw = Gx 10^ yr 



50AU 



a \ 
/im / 



Pg 



g cm 



na 



Mo 



(A9) 



where {T[a?/Csw) = l/Qsiy, where Qsw is the s tellar wind coup l ing co efficient, the ratio of 
the effective to geometric cross-sections. While iPlavchan et al.l (120051) assumed Q.s'iy = 1, 
this assumption is confirmed with modeling of dust grains in iMinato et al.l (120061 ). 



From Equations A8 and A9 it follows: 



T SW o Qrad Mg 



© 



'TPR Qsw Lq 



(AlO) 



as in 


Plavchan et al. 


Plavchan et al. 


( 


2005) 



051 ) with Qrad/ Qsw = 1 assumed. We note that Equation 3 in 
Plavchan et al.l (l2005l ) is more fundamental than Equation AlO since the solar wind mass- 



lo ss rate is vari a ble, b ut this is a good approximation. Using the X-ray mass-loss rate relation 
in 



Wood et al.l (120051 ). we can re-express AlO as: 



TSW 



Q 



rad 



TpR Qsw Lq \ L 



1.34 



0.68 



(All) 



For the Sun, Lx/Lq ~ 5 x 10^''. When tsw/'T'pr < 1, P-R drag time-scales are longer 
than stellar wind drag time-scales, and stellar wind drag is relevant. Using All and the 
approximation that 1.34=4/3, stellar wind drag dominates P-R drag when: 



Lx 
Lit 



> 3 



0.75 



> 1.13 X 10- 



m 



(A12) 



Equations All-12 are subject to the condition that Fx<8xl0^ ergs/s/cm^, which can be 
re-written as: 



< 1.26 X 10"^ 



Tq 



(A13) 



For "saturated" stars with Fx>8xl0^ ergs/s/cm^, stellar wind mass-loss can still dominate 
P-R drag, but the mass-loss rates can't be estimated accurately from the X-ray activity. 



St ellar wind drag and P-R drag are additive, and can be combined as done in lStrubbe fc Chiang 



( 120061 ): iFitzgerald et al.l (120071 ): 
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where Pcpr is defined in lStrubbe fc Chiang (120061 ): 



PcPR = Qrad + QsW 



which we can re-express as: 



Pcpr — Qrad 1 



TPR 



(A15) 



(A16) 



Drag forces predict a disk with constant vertical optical depth 
a steady state, the mass must cross an annulus is constant, 



crossmg 



constant = 27rr S^fr) Vr 



surface density). In 



(A17) 



The average radial inward drift velocity of a grain goes as Vr = r/rdrag oc r~^. Hence, Sr(r) 
must also be constant as a function of orbital radius. An infrared spectra consistent with a 
power-law follows directly from this. 

While stellar wind drag and P-R drag forces can have similar magnitudes, the same is 
not true for the stren g th of the radial or "blowout" force unless {M^/ Mq){Lq/ L^) ~ 3000 
(jPlavchan et al.l l2005l : iGustafsonI Il994r). Such high mas s-loss rates are not supported by 
recent theoretical arguments (IHolzwarth fc Jardind 120071 ). The additional factor of ~1000 
arises since the ratio of drag forces compared to the ratio of radial forces for P-R and 
stellar wind is enhanced by a factor of {c/vsw) where Vgw is the stellar wind velocity. 
Strubbe fc Chiang (l2006l ) writes the grain blowout radius from the combined effects of radi- 
ation and stellar wind as: 

3 L^PswR 



at 



where 



87r GM^cpg 

M^cVs 



PsWR = Qrad + QsW 



(A18) 



(A19) 



We estimate the minimum mass of the parent bodies - asteroid and comet proto- 
planetary analogs - by considering how qui ckly the dust mass must be replenished to main- 
tain a steady-state as in lChen et al.l (j2005al ). Grain-grain collision time-scales can be shorter 
than both P-R and stellar wind drag time-scales. The shortest time-scale dictates which 
mechanism is responsible for removing dust, but that shortest time-scale will never be longer 
than the P-R drag time-scale for the orbital dust distances typically under consideration. 
Hence, we use the P-R drag time-scale to estimate minimum total parent body mass, but 
note that the total parent body mass can be substantially larger if the grain-grain collision 
time-scale or the stellar wind drag tim e-scale are significantly shorter (jChen et al.l l2005al : 



Plavchan et al. 2005: Minato et al. 



20061 ). The minimum total parent body mass is given by: 
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MpB > ^LirU/c^ 



(A20) 



as m 



Chen &: Jural (120011 ) where t.,, is the age of the star. 



Fina lly, throughout our analys is, we assume M^, to be the average mass-loss rate. As 
noted in lAugereau fc BeustI (120061 ). flaring activity can substantially enhance the stellar 
wind drag force relative to a quiescent state. For a star that flares ~10% of the time with 
mass-loss rated enhanced by a factor of ~100 during flared states with respect to quiescent 
activity, the average mass-loss rate can be ~10 times greater than the quiescent mass-loss 
rate, and hence grain removal time-scales can be ~10 times shorter. The stellar activity 
cycle can be assumed to be much shorter than the ~Myr drag-force removal time-scale, and 
hence individual grain li fetimes will be inse nsiti ve, on average, to the flared state of a star. 
However, the methods in lWood et al.l (120051 ) and lSchmidt et al.l (120071 ) to estimate mass-loss 
rates are sensitive to the flared state of a star, and a correction factor based on the activity 
level of the star should be applied to estimate the average mass-loss rate. 



A. 2. 3. The Transition Between Grain- Grain Gollisions and Drag Forces in Disk 

Dynamical Evolution - Theory 



From Equations A14 and A6, we derive: 

0^ (l^^)iiF^j [irJ [wJ 



Tdrag 1000 / r \ 1/2+" / Q,,u \ ( Md V ( L,\ ' ( M,^^'^ 



Tcoii 7 V50AU 



We consider the dust mass estimate from the fractional infrared excess (§4.4; IChen et al. 



2005al : I.Tura et al.lll995h : 

1 r 

M, > yvr-^P.^?nin(«) (A22) 

If we set r = Dmin = Turing, and if we assume the dust is originating in a ring with relative 
width of 0.1, then Equation A21 approximates the local dust mass in the ring. Furthermore, 
Equation A14 represents the total time for a dust grain to spiral into the star. The average 
time for a dust grain to escape the ring by being dragged inwards before colliding with 
another grain is further reduced by a factor of ~ ADrmg/i^Drmg) = 0.05. If we assume 
(a) = l/xm, pg = 2.5 g cm~^, and a = 0, we get: 

Tdrag ^ . Lm/L^ ( Dring\^^'^ ( Qcoll \ / \ ^ f \ / a r,o\ 

'™g) = 777777^7^ 7777^ 7^ ^ { T^- ) (^23) 



TcoU 



For a solar-type star with a young Kuiper-belt analog at 50 AU and Mg^ = Mq such that 
Qcoll/ PcPR = 0.75, Equation A22 predicts that debris disks are drag-force dominated (and 
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collision-force dominated otherwise) when 

^<1.2x 10~^ (A24) 

Increased stellar winds increase the fractional infrared excess luminosity transition in Equa- 
tion A23 by a factor of ~ |(1 + ^{Mgw/ Mq)) for a solar-type star. If M^^ = lOM© as might 
be reasonable for young, low mass stars, the transition fractional infrared excess luminosity 
is 4x10"^ 



A. 2. 4- Planet- Dust Interaction 



When discussing the role of different grain removal mechanisms in debris disks, the 
effect of planets on the dynamics cannot be neglected. Dust trapped in orbital resonances 
with planets are invoked to explain clumps and n on-azirnuthal l y symmetric disk st r uctures 



seen i n resolved disks su ch as Vega and Eps Eri (IWyattI l2003l : iGreaves et al 



2005: Krivov 



20071 : iKrivov et al.l 120071 ). Herein we consider the dynamically cold model of IWyattI (120031 ). 
with planetesimal (and dust) eccentricities of <1%, to approximate when the presence of a 
planet is important for disk dynamical evolution relative to the processes already discussed. 



WyattI (120031 ) considers a migrating planet with migration rate hp encountering plan- 



etesimals to estimate resonance capture probabilities. Planetesi mals wil l be c aptured in 
resonances with a probability Pr given by the functional form from Iwyatti (2003): 



Pr 



H -1 



(A25) 



where do.s is the 50% capture proba bility ni i gratio n rate and 7 dictates how fast the proba- 
bility turns over as a function of dp. IWyattI (120031 ) determined that do.5 takes the form: 



•^0.5 




(A26) 



depending on the orbital reso nance (dn 5 in u nits of AU/Myr; D^; is the planet orbital radius). 
For the 2:1 orbital resonance, IWyattI (120031 ) finds through computer simulations X=5.8 and 
u=1.4; for the 3:2 resonance, X=0.37 and m=1.4; for the 4:3 resonance X=0.23 and u=1.42; 
and higher order resonances have weaker capture probabilities (bigger X). If dp < do.5, 
planetesimals will have a >50% of being caught in an orbital resonance and vice-versa. 



We relate t he planet-planetesirn al interaction analysis of IWyattl (120031 ) to planet-dust 



interaction as in 



Krivov et al.l (120071 ). First, we replace the planetesimals with dust grains. 
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Dust grains will orbit on slightly larger orbits than the Keplerian orbits of the planetesimals 
due to the effect of radiation pressure. As a result, resonant orbital semi-major axes will be 
slightly larger and grain-size dependent. The correction factor to a resonant semi-major axis 
is (1 — where /3 is the ratio of the radiation pressure force to gravity and is dependent 

on grain properties. 

Second, inward migrating dust under the action of a drag force is equivalent to an 
outward migrating planet. We can replac e with the dust inward drift velocity v^. when 
\vr\ » |dp|. In the example of Vega in IWyattI (120031 ) . a planet migration rate of dp = 
0.45AU/Myr is derived to model the resolved Vega sub-mm excess. For comparison, from 
A. 1.2, the dust migration rate in units of AU/Myr from the combined effects of P-R and 
stellar wind drag goes as: 
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(A27) 

This is a factor of ~50 larger than the inferred planet migration rate for Vega for the default 
values of grain size, density, orbital distance, etc. Hence, replacing the planet migration rate 
with the dust inward drift velocity is appropriate. 

With \vr\ » I dp I (and consequently Tcou/'Tdrag > 1), the disk would be drag- force 
dominated excepting for the presence of a planet. We expect the dust to spiral inwards 
towards the sublimation radius unless it is trapped in an orbital resonance, or scattered 
from the system in the "resonance overlap region" . 

Dust will be trapped into an orbital resonance with >50% probability if 



Vr < do.5 



1 /MpMoV M, 



D 



pi 



(A28) 



which we can re-express as: 
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If we set u=lA, we can re-write this in terms of Jovian masses (Mj): 



(A29) 
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In Equations A28 and A29, Dpj^q.q corresponds to the semi-major axis of the orb i tal re sonance 
p+q:p, which can be related to the planet semi- major axis as in iKrivov et al.l (120071 ): 



D 



p+q:q 



D 



pi 



1/3 / P + Q 

p 



2/3 



(A31) 



From Equation A29 we conclude that a proto-planet >0.05Mj at 50 AU will trap a 
majority of 1/im dust (2.5 g cm~^) around a solar type star in the 3:2 orbital resonance 
{Msw/Mq = 1 assumed). At 10 AU, the minimum planet mass is 0.14 Mj, and at 5 AU 
the minimum planet mass is 0.22 Mj. Since these planet masses are comparable to known 
extrasolar planet masses, we conclude that dust trapped in orbital resonances with Jovian 
planets should be a common phenomenon. In our own Solar System, reson ant trapping of 



a sm a ll fraction of inwa rd migrating dust has been observed for the Earth (jPermott et al. 



19941 : iReach et al.l Il995l ) , and would be observable for the Jovian planets if we were not 
embedded in our own zodiacal cloud. 



Krivov et al.l (120071 ) goes into more detail about what happens to the dust once it is 
trapped into resonance. The resulting observability is determined by the survival time of 
grains in resonance relative to the collision and drag time-scales. While in resonance, the 
dust eccentricities will be pumped up until the dust escapes resonance or collides with other 
dust. If the dust escapes resonance, it can be scattered inward and continue spiraling into 
the star, or be scattered out of the system, and this is again dependent on the planet mass. 
Additionally, the orbits of dust quickly become chaotic for planets of a few Jovian masses, 
providing an upper limit to the pl anet mass for observable r esonant trapping in addition 
to the lower limit in Equation A29 (iKuchner &: Holmaru |2003| ) . We note that the minimum 
planet mass required to capture dust is proportional to Pcpr- If th^ stellar wind mass- loss 
rate is substantially larger than solar, it affects the planet-dust interaction and disk dynamics 
overall. Stars with high mass-loss rates are less likely to have observable dust in resonance 
with planets of a similar mass compared to stars with small mass-loss rates. 

Dust in non- resonant orbits can enter the "resonance overlap region" and be scattered 
out of the system on short time-scales: 





< (^) 


^-1 


Dpi 





2/7 



(A32) 



where r and Dpi are the dust and planet semi-major axes respectively, and Mp and M^, are 
the planet and stellar masses in the same units. Again, if the planet is massive enough, the 
planet can effectively scatter a majority of the inspiraling dust out of the system, creating a 
disk with an inner hole. This is the "disk with inner hole" model is widely used to explain 
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infrared spectra with no short-wavelength "hot excess" , but longer wavelength "cold excess" . 
Such a disk would have a constant surface density truncated at Rj^ = Dpi and Rqut equal to 
the orbital radius of the colliding parent bodies producing the dust. Conversely, the presence 
of a warm excess where coUisional time-scales are long relative to drag time-scales implies 
the lack of a planet of sufficient mass interior to where dust is being produced to halt the 
inward migration. 

In summary, if Tcou/T^rag < 1, then the radial inward drift velocity of dust is effectively 
zero. Since dust grains undergo destructive collisions before they can drift into a planet, dust 
will be confined to where the parent bodies are located, such as in a ring or in resonance with 
a planet (if the planet is migrating). The presence of a planet can likely be inferred from the 
presence of clumpy disk structure. If Tcou/'Tdrag > 1, the presence of Jovian planets halts or 
substantially arrests the inward migration of dust under drag forces. If the planet is on the 
order of a Jovian mass, the dust will also likely be trapped in clumpy resonances that could 
dominate the fractional infrared excess. The presence of non-azimuthal disk structures for 
disks with L/r/L=k<10~*^ can be used to infer the presence of Jovian planets without directly 
detecting the planets themselves. Conversely, a power-law infrared excess can be used to 
exclude the presence of Jovian planets at the orbital radii of the emitting dust. 



A. 3. Spectral Type Disk Fraction Dependence 

Herein we present six factors that complicate the interpretation of the dependence of the 
70 fim disk fraction on spectral- type. Three of these factors apply particularly to our sample, 
and three are intrinsic to differences between spectral types. We consider each in turn and 
its impact on the apparent discrepancy between the 70 fim disk fractions as a function of 
spectral type. We conclude from these six factors that the paucity of bright 70/im excesses 
(F7o/F*>~5-10) around late-K and M dwarfs is likely a real phenomenon rather than due 
to sensitivity or systematic bias. 

First, the average age of the T^,<5000 K sub-sample is 51 Myr (median 12.25 Myr), 
compared to the average age of 380 Myr (median 340 Myr) for the 5000 K < T^, < 6000 
K sub-sample. Second, the average fiux ratio upper-limit of the T=i,<5000 K sub-sample is 
F7o/F,=18.5 (median 6.1), and is F7o/F,=2.1 (median 1.8) for the 5000 K < T, < 6000 K 
sub-sample. The first two factors cancel to first order, when accounting for the third factor 
- the assumed ag^e evol ution of the fractional infrared excess. We assume that L//j/L* oc t;^^ 
(e.g.. 



Rieke et al 



l2005l ). While we are ~3-9 times less sensitive on average to 70 /im excesses 
around the cooler stars, a steady-state collisional evolution model predicts that they should 
be ~10 times brighter at younger ages. 
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Note that the first t wo factors do not can cel if we assume a different time-dependence for 
the evolution of L/^j/L^,. ICurrie et al.l (l2008bl ) posits that debris disks around solar type stars 
undergo a "rise and fall" of debris disk brightness, peaking at ages of ~10-15 Myr. Low mass 
K and M dwarfs could undergo a similar but delay ed "rise and fall" of d e bris disk brightness, 
peak ing at ages of 30-50 Myr (e.g. NGC 2547, iForbrich et al.l l2008l : iKenyon fc Bromley 
20081 ) . Under such an alternate scenario, for the median age of the T^,<5000 K sub-sample, 
we would be catching the stars during the rise and before the peak in emission from debris 
disks. 

Fourth, we have a slightly higher fraction of excesses at 24 fim for stars with T* < 5000 
K (18.5^1%) than for stars with 5000 K < T, < 6000 K (14.3|°%). This result implies 
that small 24//m excesses (e.g. ~10-15% flux density excess) for young T^, < 5000 K stars 
are as common as they are around hotter stars, whereas larger fractional infrared excesses - 
e.g. >100% as probed with our 70 /im observations - remain rare around K and M dwarfs 
relative to higher mass stars. The lack of large fractional infrared excesses around K and M 
dwarfs still necessitates an explanation. 

Fifth, at 70 the emitting dust orbits at smaller radii for cooler stars. At small 
orbital radii (~3-10 AU, Table 4), we are equally sensitive to dust for all stars in our 
sample. However, the hotter stars with strong 70 fim excesses have dust at orbital radii that 
we cannot probe with the 70 /xm observations of the cooler stars. The fifth factor implies 
that the lack of low mass debris disks we observe at 70 /im could depend on disk geometry. If 
young Kuiper belt analogs (that produce the dust we are observing) form at radii that scale 
with the stellar luminosity, then we are indeed observing fewer massive, bright debris disks 
around the cooler K and M dwarfs compared to G dwarfs at comparable fractional infrared 
excess (e.g., F7o/F^, > 5-10). Theoretical modeling indicates that the orbital distance for 
the formation of Jovian p l anets at the "snow line" h as a weak dependence on stellar mass 
(IKennedy fc Kenyonll2008l : Kennedy et al.ll2007l . l2006l ). This implies that young Kuiper Belt 
analogs may form at smaller radii for lower mass stars, but this is not conclusive. 

If instead the young Kuiper Belt analogs form at comparable orbital radii independent 
of spectral type, our observations are not sensitive t o detect these cold debris disks around 
late-K and M dwarfs. Sub-millimeter observations by lLestrade et al.l (120061 ) : iLiu et al.l (120041 ) 
have revealed three cold debris disks around M dwarfs, albeit at low S/N for GJ 842.2. The 
second sub-mm excess source, GJ 182, is not detected in excess in our sample, although an 
excess 5 times the photosphere at 70 /im would not have been detected. The third source, 
AU Mic, is detected in excess of ~9 times the photosphere at 70/im in our sample. These 
initial sub-mm observations of low-mass stars imply cold sub-mm disks around late-K and M 
dwarfs could be as common as 70 /xm excesses around older G dwarfs (~15%, iBryden et al. 



20061 ). Observing more M dwarfs at sub-mm wavelengths will further address this possibility. 



The dust removal mechanisms and time-scales are the sixth factor. The dust lifetimes 
around K- and M-type stars due to P-R drag and collision time-scales are (sometimes signifi- 
cantly) longer and the dust blowout radii are smaller, scaling with either the stellar luminosity 
or the stellar mass. For example, an MO dwarf has a P-R drag time-scale ~10 times longer 
than for a G2 star, a collision time-scale 40% longer (for the same disk mass), and a blowout 
radius for dust grains ~10 times smaller (these follow from §A.l). A coUisional-cascade 
steady-state predicts the number density of grains scales as a^'^'^ where a is the grain size 
and the minimum grain size is set by the blowout radius. The smaller blowout radius and 
longer grain lifetimes for K and M dwarfs imply larger dust abundances and more smaller 
grains relative to G-type stars for a similar amount of colliding parent bodies. Such an excess 
of dust should be more easily detectable, and this is not observed at 70 fim. A mechanism 
like stellar wind drag is necessary to account for this paucity of bright 70 fim debris disks 
around M dwarfs. 
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Table 1. Stellar Properties 



Name 


HD / 


Spectral 




Assoc ROSAT 


Lithium 


log(R^^)'^ Stellar Rotation References'^ 




Name 2 


Type'' 


(pc) 


-iation'' (cnts/sec) 


Abundance 


vsmi Period 












log(N(Li)) 


(km/s) (days) 












(H=12) 







1835 " 


GO 


20 


HSC 


2.89E-01 


2.58±0.10 


-4.46 


3.6 


7.81 


1,1,2,1 




7590 


GO 


24 




1.91E-01 


2.75±0.10 


-4.46 


6.4 


5.67 


1 


BB Scl'^'S 


9770 


G5 


24 




2.59E+00 






2 




3 




10008 


G5 


24 


LA 


1.44E-01 


1.96±0.12 


-4.38 


2.9 




1 




11131 


GO 


23 




3.86E-01 


2.42±0.10 


-4.47 


3.4 


8.92 


1,1,4,1 




15089f'9 


A5p 


43 




2.28E-01 






49 


1.74 


5,6,7 




16287'^ 


KO 


24 






0.55 


-4.32 


4.9 




8 




16555 


A5 


44 




3.20E-03 






315 




9,7 




17240^ 


F2 


47 










45 




10,11,12 


DO Eri 


24712 


FO 


49 










18 


12.46 


5,6,11,12 


v833 Tau« 


GJ 171.2 


K3' 


18 




2.69E+00 




-4.057 






13 




29697 


K2 


13 


UMG 


2.36E+00 


0.754 


-4.1 


9.5 


4 


14,13,15,15 


Gl 182 




MQj 


27 


LA 


6.51E-01 


1.8 




14 


1.9 


16,3,17 




36435 


GO 


20 




1.42E-01 


1.60±0.15 


-4.44 


4.5 




1,1,18 


AB Dor"' 


36705 


GO 


15 


ABDMG 


7.22E+00 


3.1 




53 


0.51 


16,3,17 




41824^ 


G5 


30 




1.90E+00 




-4.17 


8 


3.3 


19,3,20 




41593 


KO 


15 


UMG 


2.57E-01 


1.22±0.17 


-4.36 


5.0 


7.97 


1 


GJ 3400Af 


48189 


GO 


22 


ABDMG 


2.33E+00 


3.3 


-4.29 


17.6 




16,19,3 




51849'' 


KO 


22 




5.03E-01 














52698 


KO 


15 




1.99E-01 




-4.59 


3.5 




1 




59967 


G5 


22 




3.20E-01 


2.50±0.10 


-4.44 


3.4 




1 




63433 


GO 


22 


UMG 


2.24E-01 


2.33±0.10 


-4.34 


6.1 


6.46 


1 




72760 


G5 


22 




9.65E-02 


0.57±0.45 


-4.39 


4.2 




1 




73350S 


GO 


24 


HSC 


1.30E-01 


2.09±0.10 


-4.49 


3.8 


6.14 


1 




74576 


KO 


11 




4.60E-01 




-4.31 


3.9 




19,3 




92139' 


F2 


27 




2.03E-01 






15 




9,11,12 




92945 


KO 


22 


LA 


1.18E-01 


2.41 


-4.39 


5.1 




8,21 




97334S 


GO 


22 


LA 


3.20E-01 


2.64±0.09 


-4.41 


5.6 


8.25 


1 


TWA 14'^ 




MO 


55 


TW Hya 


1.07E-01 








0.63 


22 


TWA 12l^ 




M2 


55 


TW Hya 


4.29E-01 






15 


3.303 


3,23 


TWA ISA"^ 




Ml 


55 


TW Hya 


1.16E-01 






12.3 


5.56 


3,22 


TWA 13B'' 




Ml 


55 


TW Hya 


1.16E-01 








5.35 


22 



Table 1 — Continued 



Name 


HD / 


Spectral 




Assoc 


ROSAT 


Lithium 


log(R^^)'^ Stellar Rotation References'^ 




Name 2 


Type^ 


(pc) 


-iation'' 


(cnts/sec) 


Abundance 


vsini Period 














log(N(Li)) 


(km/s) (days) 














(H=12) 







103928° 


FO 


47 




5.00E-03 






94 


0.70 


5,24,11,12 




105963S 


KO 


27 




1.50E-01 


2.06 


-4.22 


6 


7.44 


8,8,8,20 




109011 


KO 


24 


UMG 


6.00E-02 


0.81±0.15 


-4.36 


5.5 


8.81 


1 


TWA ISB'' 




M2 


55 


TW Hya 


1.16E-01 








0.72 


22 


TWA ISA'' 




Ml. 5 


55 


TW Hya 


1.16E-01 








0.65 


22 


TWA le"" 


112429 


Ml. 5 

FO 


55 

29 


TW Hya 


1.07E-01 






130 




9,11,12 




113449 


KO 


22 




1.70E-01 


2.08±0.12 


-4.34 


5.8 


6.54 


1,13,1,1 


TWA 17^ 




K5 


55 


TW Hya 


7.21E-02 








0.69 


22 


TWA 18'^ 




MO. 5 


55 


TW Hya 


6.22E-02 








1.11 


22 




116956 


G5 


22 


LA 


2.01E-01 


1.42±0.12 


-4.447 


5.6 


7.8 


1,13,1,1 




123998 


A2p 


43 










15 




9,11,12 




125158 


A3 


48 














11,12 




128987 


G5 


24 




9.81E-02 


1.78±0.11 




3.2 


9.35 


1 




128400 


GO 


20 


UMG 


1.40E-01 




-4.56 


4.4 




1,3 




135599 


KO 


16 


UMG 


9.93E-02 


0.34±0.28 




4.6 


5.97 


1 




141272 


G5 


21 


LA 


1.87E-01 


0.48±0.27 


-4.452 


4.0 


14.01 


1,13,1,1 




144197 


A3p 


38 










7 




25,11,12 




148367^ 


A2 


37 




2.88E-01 






18 




25,11,12 




165185 


G5 


17 


UMG 


5.85E-01 


2.65±0.10 


-4.49 


7.2 




1 


RE1816+5411 


EY Dra 


Ml. 5 


30 


LA 


2.92E-01 






61 


0.459 


26,17 




175742 


KO 


21 


UMG 


1.88E+00 


1.39 


-3.90 


14 


2.907 


8,8,8,23 




177724 


AO 


26 




1.41E-01 






317 




5,11,12 




180161 


KO 


20 




1.55E-01 


0.86±0.14 


-4.44 


3.1 


9.7 


1 




186219 


A3 


42 










125 




9,7 


AT Micf 


196982 


M4.5j 


10 


BPMG 


3.91E+00 






10.1 




3 


eta Ind 


197157 


FO 


24 














11,12 


AU Mic<= 


197481 


MOJ 


10 


BPMG 


5.95E+00 






9.3 


4.9 


3,17 


BO Mic 


197890 


KO 


44 


LA 


6.11E+00 






120 


0.38 


27,28 




358623™ 


K6 


24 


BPMG 


2.36E-01 






15.6 


3.41 


3,23 


SAO 145139^ 




K5j 


26 




6.03E-01 














202730° 


A5 


30 




4.19E-01 






135 




25,11,12 
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Name 


HD / 


Spectral 




Assoc 


ROSAT 


Lithium 




Stellar Rotation 


References'^ 




Name 2 


Type'' 


(pc) 


-iation'' 


(cnts/sec) 


Abundance 




vsini 


Period 
















log{N{Li)) 




(km/s) 


(days) 
















{H=12) 












203244 


G5 


20 




1.90E-01 


1.64±0.15 


-4.39 


4 




1,1,3 


LO Peg 


HIP 106231 


K3' 


25 


ABDMG 


9.23E-01 


1.624 


-3.906 




0.42 


14,13,23 


Gl 859ABf 


212697 


GO 


16 




8.47E-01 


2.9 




9.7 




29,18 


HKAqr 


Gl 890 


MOJ 


22 


CMG 


4.27E-01 








0.431 


17 




218738S 


G5 


25 




1.93E+00 


2.13 


-3.90 


7 




8 


Gl 907. if 




K7j 


27 




4.04E-01 






3.4 




3 



'^Spectral type from the HD Catalog unless otherwise noted dNesterov et al.lll995l ). Trigonometric parallax distance from Hipparcos 
unless otherwise noted. 

''TW Hya = ~ 8 Myr T W Hya association l lKastner et al.l Il997l : IWebb et all Il999l): BPMG = ~1 2 Myr /3 Pic moviii g 
group llZuckerman et al.ll200lh : ABDMG = ~50-100 Myr A B Dor moving group llZuckerman et al.l [iool iLuhman et"al] l2005h : 
LA = ~20-100 Myr Loca l Association llMontes et al.l l2001ah : CMG = ~200 Myr Castor moving group iMontes et alTl2001al : 



LA = ~20-100 Myr Loca l Association llMontes et al.l I2001af) : CMG = ~200 Myr Castor moving group llMontes et al.l I2001al : 

iBarrado fc Navascueslll998h: UMG = ~30 0-500 Myr Ursa major group l lKine et al.ll2003l : ISoderblom\fe Mavoj|l993al lbl): and HSC ^ 
= Hyades supercluster lMontes et al.ll2001al ) 

'^Calcium II R^^ = log{h jj K /^bol) corrected for photospheric contribution. 

"^References listed in order of Association, Lithium Abundance, log(R^^), Stellar Rotation, and young A star sources where appro- 
priate. 

"Visual companion with separation >5" not detected with Spitzer. The visual companion is of unknown age/distance, and is 
excluded from our sample. Projected separation is greater than 100 AU. 

'All photometry in this paper includes both A and B components (separation <5"), including co-adding photometry if the binary 
is resolved in some (but not all) bandpasses. Projected separation is less than 100 AU. 

^Visual companion with separation >5" detected with Spitzer at 24 fim and not at 70 /im. The visual companion is of unknown 
age/distance, and is determined to be photospheric at 24 /im from the Ks-[24] color. Projected separation is greater than 100 AU. We 
exclude this companion from our sample. 

•"Einstein X-ray detection of 3.7x10"^'' ergs/cm^ /s l |johnson|ll986l '). 

iGrav et al.l ||2003|) 

■ iHawlev et al] lll996h 



Distance set to media n distance of TWA 1-19 members with known trigonometric parallax as in iILow et al.ll2005h . We note that 
iLawson &: Crause ] ll2005l) present evidence for a median distance of ~90pc. Due to uncertaint ies in age and distance, we do not use 
photometric distances on a star by star basis. Spectral types from lSarrado &c NavascuesI ll200d) . 

'Distance estimated from absolute K magnitude dis tance modulus for a 115 Myr Ml. 5 dwarf. Age inferred f rom p ossible Pleiades 
measurement and hx /L* implied age of <300 Myr l lMessina et al.ll2003h . Spectral type from [jeffries et al.l l ll994) . The distance 
derived from the Ks magnitude is not consistent with Pleiades membership; hence the age, distance, and estimated radius for this 
source are not well-constrained. 

■"Distance estimated from absolute K magnitude distan ce modulus for a 12 Myr K7 dwarf. Age inferred from membership in the /3 
Pic moving group. Spectral type from Flbrres et al] | |2006|) . 



References. — 1 - ICaidos ct al.' bOOOl). 2 -IValenti fc Fischer] Jiooj) . 3 -iTorres et al.l ll200d). 4 -ISoderblom fc Mavoj (^1993a^, 5 - 
iKudrvavtsev fc Roma nvuk (2003), 6 - iRover et al. I 20021) 7 - Stauffei l|priy.comni.h . 8 - fstrassmeier et al.l 1 I2OOOI) . 9 - Glebocki et al.] 

' Sond boodi ). 13 - iGrav et al.l l20oi) . 14 - iFischei (1998|), 15 

Messina et al.l ||200;), 18 - IPace et al. ||200;), 19 -iHenrv et ah J 19991. 20 
Lawson fc Crause 

iBromaee et al.l lll992h . 28 - ICutispoto et al.l lll997^ . 29 - 



ll2000h. 10 -iDanziger fc FabeJ (11973). 11 - ISong et al.l 



iMontes et al 



ll2001bl). 16 - iFavata et al.l ( Il998i) . 17 



-iKazarovets et al.ll2006l). 21 - 
(l2002h. 25 -lErspamer fc Nort 
IPallavicini et al.l \l98% 



Grav et al.l ( l2006h. 22 
||2003|) . 26 - Ijefltries et al.l (Il994h . 27 
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Table 2. MIPS 24 /^m and 70 /im Fluxes (Not Color- Corrected) 





Measured'' 


Predicted 




Measured'' Predicted 










MIPS 


Photosphere'' 




MIPS Photosphere*' 






New 




F^{24 fim) 


F^(24 fim) 




F^(70 fim) F^{70 /im) 






Detection of 


Name/HD 


(mJy) 


(mJy) 


X24 


(mJy) (mJy) 


X70 


References'^ 


Excess? 


Stars with MIPS- 24 and MIPS- 70 excesses'^ 


10008 


41.6±0.5 


36 


12.1 


30±1« 


4.1 


19.8 


1 


Y 


73350 


66.2±1.3 


58 


6.22 


123±4 


6.6 


32.3 




Y 


112429 


132.8±1.6 


121 


7.11 


58±4* 


14 


10.3 


1 


N 


135599 


82.5±0.9 


74 


9.15 


106±10 


8.5 


9.4 




Y 


Stars with MIPS- 70 excesses'* 


7590 


62.7±1.0 


61 


2.07 


233±16 


6.9 


13.9 




Y 


59967 


70.5±2.5 


65 


2.14 


36±4 


7.4 


7.5 




Y 


92945 


39.5±1.0 


39 


0.57 


305±15 


4.4 


19.4 


1 


N 


AU Mic 


155.2±3.2 


150 


1.65 


223±26 


17 


8.0 


1,2 


N 


Stars with MIPS-24 excesses'* 


GJ 3400A 


120.1±1.4 


109 


8.25 


16±5S 


12 


0.7 


1 


Y 


123998 


136.6±2.3 


123 


5.7 


<20 


14 






Y 


175742 


46.9±0.6 


41 


10.8 


<13 


4.7 




1 


Y 


AT Mic 


131.6±2.lf 


114 


8.21 


22±5* 


13 


1.8 


1,2 


Y 


BO Mic 


16.4±0.4 


14 


5.66 


<12 


1.6 




1 


Y 


358623 


14.2±0.2 


13 


7.28 


<9 


1.5 




1,2 


Y 


Gl 907.1 


26.7±0.4 


24 


5.44 


<13 


2.8 




1 


Y 


Stars with no MIPS excesses'* 


1835 


86.0±1.0 


81 


4.95 


18±4S 


9.3 


2.3 


1 




BB Scl 


87.6±1.4 


97 


-6.89 


<13 


11 




1 




11131 


64.2±2.3 


62 


0.83 


<16 


7.1 




1 




15089 


142.1±8.6 


143 


-0.08 


28±3'' 


16 


3.7 






16287 


30.8±0.8 


30 


0.89 


<13 


3.4 








16555 


108.5±1.7 


111 


-1.24 


<12 


13 








17240 


56.2±0.7 


54 


3.07 


31±6* 


6.2 


3.9 






DO Eri 


56.6±1.0 


56 


0.53 


<12 


6.4 








v833 Tau 


63.5±2.5 


62 


0.77 


<28 


7.0 








29697 


69.7±0.7 


68 


2.53 


<31 


7.7 




1 




Gl 182 


28.7±0.7 


27 


2.06 


<19 


3.1 




1 




36435 


61.0±0.7 


59 


2.24 


<12 


6.8 








AB Dor 


106.0±8.8 


95 


1.22 


<47 


11 




1 




41824 


85.4±1.5 


84 


0.74 


<11 


9.6 








41593 


84.6±1.4 


84 


0.34 


<22 


9.6 








51849 


27.9±0.4 


26 


4.67 


<10 


3.0 




1 




52698 


102.1±2.0 


100 


1.15 


<13 


11 








63433 


58.6±0.6 


56 


3.64 


<13 


6.4 








72760 


48.8±0.8 


48 


0.61 


<12 


5.5 








74576 


130.2±1.7 


129 


0.71 


<16 


15 




1 




92139 


419.5±3.5 


403 


4.57 


59±18 


46 


0.7 


1 




97334 


73.3±1.3 


74 


-0.66 


<15 


8.5 








TWA 14 


4.21±0.20 


3.8 


2.24 


<13 


0.43 




3 
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Table 2 — Continued 





Measured^ 


Predicted 




Measured^ 


Predicted 










MIPS 


Photosphere'' 




MIPS 


Photosphere'' 






New 




F^(24 fim) 


F^(24 /xm) 




F^(70 /im) 


F^(70 fim) 






Detection of 


Name/HD 


(mJy) 


(mJy) 


X24 


(mJy) 


(mJy) 




References'^ 


Excess? 


TWA 12 


6.04±0.33 


5.7 


1.17 


<8 


0.65 




3 




TWA 13A 


8.67±3.44 


9.1 


-0.11 


<38 


1.0 




3 




TWA 13B 


8.67±3.44 


9.0 


-0.09 


<38 


1.0 




3 




103928 


42.0±1.0 


41 


0.83 


<11 


4.7 








105963 


36.1±11.8 


36 


-0.03 


<10 


4.2 




1 




109011 


39.7±0.8 


40 


-0.13 


<11 


4.6 








TWA 15B 


1.60±0.62 


1.4 


0.31 


<76 


0.16 




3 




TWA 15A 


1.44±0.62 


1.3 


0.27 


<76 


0.15 




3 




TWA 16 


5.87±0.14 


5.5 


2.91 


<17 


0.62 




3 




113449 


48.5±0.8 


45 


4.34 


<15 


5.2 








TWA 17 


2.20±0.80 


2.1 


0.17 


<18 


0.24 




3 




TWA 18 


2.78±0.09 


2.7 


0.64 


<9 


0.31 




3 




116956 


50.5±1.7 


49 


0.9 


14±2f 


5.6 


3.6 






125158 


104.3±9.8 


97 


0.77 


<429 


11 








128987 


44.3±0.5 


44 


1 


<22 


5.0 








128400 


67.0±1.2 


67 


0.15 


<11 


7.6 








141272 


43.3±0.8 


45 


-2.13 


<14 


5.1 




1 




144197 


131.8±1.4 


140 


-6.04 


<26 


16 








148367 


159.0±2.2 


154 


2.25 


<23 


18 








165185 


115.1±3.8 


116 


-0.35 


<16 


13 








RE1816+541 


7.52±0.13 


7.4 


0.72 


<13 


0.85 




1 




177724 


498.2±8.5 


485 


1.56 


50±14f 


55 


-0.4 


1 




180161 


59.2±0.6 


59 


-0.14 


<11 


6.8 




1 




186219 


90.2±0.9 


86 


4.38 


23±6 


9.8 


2.1 


1 




eta Ind 


201.3±3.1 


212 


-3.37 


31±6 


24 


1.2 


1 




SAO 145139 


21.1±0.9 


21 


0.37 


<15 


2.4 








202730 


157.3±17.1 


157 


0.02 


27±4S 


18 


2.2 


1,4 




203244 


61.2±1.1 


59 


2.47 


<14 


6.7 








LO Peg 


23.6±0.4 


21 


4.97 


<8 


2.4 




1 




Gl 859A 


165.3±2.9 


163 


0.7 


<18 


19 




1 




HK Aqr 


13.3±0.4 


13 


0.29 


<16 


1.5 




1 




218738 


44.1±10.4 


39 


0.5 


<8 


4.4 




1 





^Flux from DAT pipeline PSF fitting unless otherwise noted. Detection confirmed with MOPEX pipeline PRF fitting 
unless otherwise noted. 

''From Ks-[24] empirical model derived in iGautier et all l|2008l ) and stellar effective temperature. Stellar effective 
temperatures derived from PHOENIX NextGen model photospheres fit to optical and near-IR photometry. 

'^Previous publications of Spitzer MIPS photometry for these sources. 

''A 70 /xm excess is defined as X7o>5. A 24 /xm excess is defined as X24>5. X24,70 = (^24, 70 — F',)/a24,jo where the 
uncertainty 0-24,70 does not include the uncertainty in the model photosphere flux density, but does include the uncertainty 
in the detection S/N and the flux density calibration uncertainty. 

°Centroided position of 70 /im detection is offset 3. "7 in right ascension and 2. "9 in declination from position expected 
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from 24 fim detection. 

'DAT aperture corrected flux density. Detection confirmed with MOPEX PRF fitting. 

spiux from MOPEX pipeline PRF fitting. Detection not confirmed with DAT pipehno PSF fitting. 

'^DAT pipeline detection not confirmed with MOPEX pipeline PRF fitting. 

'Centroided position of 70 fira detection is offset 7. "4 in right ascension and -l."5 in declination from position expected 
from 24 fim detection. 

References. — l 4Chen et ahl ll2005ah . 2 4Rebull et ahl bOO^ . 3- ILow et al.l l l2005lV 4 4Rreke et al.l boodt) 



Table 3. Optical and Near-IR Photo-electric Photometry and Derived Stellar Properties 



Name/HD 




T*^ 




U 


B 


V 


Rc 


I 




J 


H 




References 
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''Radii are fit to the nearest 0.5% and the uncertainty is dominated by the uncertainty in the trigonometric parallax. For unresolved binaries listed 
in Table 1, the derived radius, effective temperature, and luminosity listed correspond to the equivalent single star (e.g., h^,=hj^+hg; R^^v^ Ra if 
Tt~Tyi~Ta and R^^R^). For sources with non-trigonometric distances in Table 1, the derived radius and luminosity are dependent on the distance 
assumed. 

''Not included in fit to model photospheres. 
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Table 4. Derived Stellar and Thin Dust Disk Model Dust Properties 



Name/HD 








log(Ljfl/L,)b 


Md MpB 


Iog(Lx/L*) 


Man, 


Age 


Age 


M, 




(mjy) 


(K) 


(AU) 




(Minoon) (Mnioon) 






(Myr) 


Method 


(Mo) 


Stars with MIPS-24 and MIPS-70 excesses 


10008 


33±1 


<78 


>8.7 


-4.09 


2.4E-05 >0.11 


-4.44 




275±25 


1,3 


0.94 


73350 


136±4 


<63 


>19 


-3.84 


3.7E-04 >0.6 


-4.80 




513±136 


1,2,3 


1.09 


112429 


63±5 


<86 


>24 


-4.69 


3.2E-04 >0.2 


<-5.62 




250±200 


5 


1.56 


135599 


117±12 


<62 


>12 


-3.88 


6.6E-05 >0.05 


-5.11 


14 


214±17 


3,4 


0.92 


Stars with MIPS-70 excesses 


7590 


259±18 


<40 


>51 


-3.57 


5.2E-03 >1.1 


-4.69 




462±35 


1,2,3 


1.15 


59967 


39±4 


<71 


>14 


-4.44 


4.6E-05 >0.12 


-4.46 




353±68 


1,2 


1.09 


92945 


341±17 


<46 


>23 


-3.12 


1.2E-03 >1.0 


-4.47 




294±23 


1,2 


0.93 


AU Mic 


247±28 


<59 


>6.8 


-3.36 


2.0E-03 >0.004 


-2.82 




12±2 


4 


0.31 


Stars with MIPS-24 excesses 


GJ 3400A 


16±5 


<200 


>2.2 


-4.86 


7.6E-07 >0.006 


-3.76 




75±25 


4 


1.03 


123998 




<100 


>31 


-5.22 


3.9E-04 >0.35 


<-5.74 




250±200 


5 


1.77 


175742 




<110 


>3.5 


-4.46,<-4.27'"^ 


1.3E-06 >0.002 


-3.18 




51±3 


3 


0.77 


AT Mic 


23±5 


<120 


>1.0 


-4.24 


3.0E-07 >0.0003 


-2.82 




12±2 


4 


0.15 


BO Mic 




<81 


>7.8 


-4.43,<-3.90'^ 


8.5E-06 >0.003 


-2.16 




25.5±8.5 


7 


0.90 


358623 




<88 


>2.7 


-4.32,<-3.70<^ 


2.8E-07 >0.0003 


-3.36 




12±2 


4 


0.67 


Gl 907.1 




<93 


>3.9 


-4.48,<-3.91<^ 


1.3E-06 >0.001 


-3.47 




12.5it5.5 


7 


0.60 


Stars with no MIPS excesses 


1835 


19±4 


<86 


>11 


<-4.97 


<8.5E-6 


-4.60 




466±55 


1,2,3 


1.09 


BB Scl 




<260 


>1.1 


<-4.69 


<1.4E-7 


-3.45 




180±120 


1,7 


1.17 


11131 




<91 


>9.2 


<-4.83 


<8.8E-6 


-4.39 




506±255 


1,2,3 


1.07 


15089 


30±3 


<100 


>38 


<-5.47 


<5.1E-4 


-5.42 




200±150 


5 


1.87 


16287 




<76 


>7.9 


<-4.41 


<7.5E-6 


-4.71 


41 


340±252 


2 


0.88 


16555 




<270 


>3.7 


<-5.21 


<5.3E-6 


-6.87 




200±150 


5 


1.56 


17240 


34±7 


<62 


>48 


<-4.68 


<1.5E-3 


<-5.22 




1100±550 


6 


1.37 


DO Eri 




<98 


>22 


<-5.21 


<l.lE-4 


<-5.31 




160±110 


5 


1.58 


v833 Tau 




<73 


>7.4 


<-4.20 


<8.6E-6 


-3.09 




15±3 


3 


0.81 


29697 




<77 


>5.1 


<-4.15 


<3.0E-6 


-3.21 




49±37 


2,3 


0.75 


Gl 182 




<69 


>6.5 


<-3.78 


<1.9E-5 


-3.18 




9±4 


7 


0.45 


36435 




<100 


>5.6 


<-4.87 


<1.9E-6 


-4.79 




504±145 


1,2 


0.98 


AB Dor 




<80 


>8.0 


<-4.31 


<1.2E-5 


-3.04 




75±25 


4 


0.91 


41824 




<300 


>1.0 


<-4.71 


<2.0E-7 


-3.61 




109±9 


3 


1.38 


41593 




<90 


>6.5 


<-4.72 


<3.1E-6 


-4.59 




326±93 


1,2,3 


0.93 



Table 4 — Continued 



IN ame / nu 


r 70 
(mJy) 




/ ATT^ 




NLO 
(,iVimoon ) 


IvlpB 

/A/r ^ 
(Mmoon ) 


iOg{ljX / i-i* ) 


A/f 


Age 
(Myrj 


Age 
Method 




51849 




<94 


>3.6 


<-4.16 


<1.6E-6 




-3.48 




65±25 


7 


0.71 


52698 




<280 


>1.0 


<-4.74 


<3.1E-8 




-4.73 


48 


809±379 


1,2 


0.89 


63433 




<30 


>73 


<-4.75 


<5.1E-4 




-4.55 




299±96 


1,2,3 


1.03 


72760 




<90 


>6.8 


<-4.70 


<3.9E-6 




-4.74 


57 


415±194 


1,2 


0.94 


74576 




<300 


>1.0 


<-3.75 


<9.6E-10 




-4.43 




215±94 


1,3 


0.84 


92139 


62±19 


<120 


>21 


<-5.50 


<l.lE-4 




-5.76 




100±50 


5 


1.58 


97334 




<100 


>8.0 


<-4.96 


<5.3E-6 




-4.53 




451±199 


1,2,3 


1.12 


TWA 14 




<50 


>8.7 


<-3.01 


<1.7E-4 




-3.00 




8±2 


4 


0.27 


TWA 12 




<60 


>7.5 


<-3.41 


<7.3E-5 




-2.57 




8±2 


4 


0.27 


TWA 13A 




<50 


>14 


<-2.97 


<l.lE-3 




-3.26 




8±2 


4 


0.32 


TWA 13B 




<50 


>15 


<-2.99 


<1.0E-3 




-3.29 




8±2 


4 


0.38 


103928 




<67 


>38 


<-5.18 


<2.4E-4 




-6.23 




440±390 


5 


1.53 


105963 




<90 


>6.7 


<-4.63 


<4.7E-6 




-4.39 




266±11 


1,3 


0.87 


109011 




<87 


>6.4 


<-4.54 


<4.1E-6 




-4.65 


50 


348±120 


3 


0.84 


TWA 15B 




<37 


>9.7 


<-1.71 


<6.4E-4 




-2.41 




8±2 


4 


0.66 


TWA 15A 




<35 


>10 


<-1.63 


<1.4E-3 




-2.37 




8±2 


4 


0.38 


TWA 16 




<53 


>9.4 


<-3.07 


<1.8E-4 




-3.17 




8±2 


4 


0.38 


113449 




<76 


>8.6 


<-4.51 


<8.4E-6 




-4.26 




212±33 


2,3 


0.86 


TWA 17 




<51 


>7.1 


<-2.74 


<5.3E-5 




-3.02 




8±2 


4 


0.81 


TWA 18 




<52 


>6.8 


<-3.01 


<5.3E-5 




-2.98 




8±2 


4 


0.38 


116956 


15±3 


<86 


>7.7 


<-4.64 


<5.6E-6 




-4.34 




334±96 


1,2,3 


0.98 


125158 




<53 


>110 


<-3.87 


<1.0E-1 




<-5.61 




125±75 


5 


1.70 


128987 




<73 


>11 


<-4.44 


<2.3E-5 




-4.85 


54 


623±74 


1,3 


1.01 


128400 




<130 


>4.1 


<-5.03 


<8.8E-7 




-4.86 


62 


795±143 


1,2 


1.03 


141272 




<83 


>7.5 


<-4.59 


<5.4E-6 




-4.40 




567±382 


1,2,3 


0.93 


144197 




<110 


>26 


<-5.37 


<1.9E-4 




<-5.81 




62.5±12.5 


5 


1.75 


148367 




<120 


>21 


<-5.58 


<8.2E-5 




-5.22 




290±240 


5 


1.77 


165185 




<160 


>3.2 


<-5.14 


<5.6E-7 




-4.46 




437±186 


1,2 


1.13 


RE1816+541 




<57 


>5.1 


<-3.30 


<2.0E-5 




-2.86 




115±100 


1,4 


0.23 


177724 


52±14 


<230 


>9.4 


<-5.75 


<2.1E-5 




-6.26 




100±50 


5 


2.30 


180161 




<110 


>4.7 


<-4.92 


<1.2E-6 




-4.68 




468±61 


1,2,3 


0.96 


186219 


25±7 


<89 


>30 


<-5.11 


<3.1E-4 




<-5.55 




200±150 


5 


1.68 


eta Ind 


32±6 


<140 


>11 


<-5.41 


<1.6E-5 




<-5.91 




250±200 


5 


1.62 


SAO 145139 




<64 


>8.1 


<-3.98 


<1.3E-5 




-3.25 




35±7 


7 


0.79 



Table 4 — Continued 



Name/HD 


F70" 






log{L,fl/L,)'' 




MpB 


\og{Lx/U) 




Age 


Age 


Mh, 




(mJy) 


(K) 


(AU) 




(Mmoon) 


(Mmoon) 






(Myr) 


Method 


(Mq) 


202730 


30±4 


<59 


>76 


<-5.77 


<4.6E-4 




-5.11 




150±100 


5 


1.87 


203244 




<53 


>21 


<-4.99 


<2.2E-5 




-4.58 




334±79 


1,2 


0.99 


LO Peg 




<96 


>3.6 


<-4.15 


<1.8E-6 




-3.13 




75±25 


4 


0.71 


Gl 859A 




<270 


>1.3 


<-4.96 


<2.0E-7 




-4.48 




314±63 


1 


0.99 


HK Aqr 




<61 


>4.5 


<-3.48 


<8.2E-6 




-3.02 




25.5±12.5 


7 


0.29 


218738 




<120 


>3.5 


<-4.61 


<1.3E-6 




-3.22 




40±8 


7 


0.85 



"Color-corrected to the dust temperature. 

^Unless otherwise noted, the fractional infrared luminosity is determined from a single temperature blackbody fit to the 24 and 70 ^m 
photometry (or upper limit). 70 fj,m photometry are color-corrected to the maximum blackbody temperature as described in §4.4. 

^^Mq ~ 2 X lO^'^-^ Mq/jv 

'^The first fractional infrared excess uses the approximation L/^ ~ uFi, at 24 /im after subtracting off the expected stellar flux density 
and ignoring the 70 /im flux density upper-limit. The fractional infrared excess upper limit is derived from the 70 fira upper limit using 
the single-temperature blackbody fit as described in table note (b). 



References. — 1 - X-ray — age correlation (§4.3, iMamaiek &: HiU enbran^ [200 ^, 2 — Calcium H&K emission-age correlation 
llMamaiek fc HillenbraiidlbOOSl. §4.3 , ) , 3 - Rotation— age correlation |Mama iek fc Hillenbrand! bOOSi , §4.3, ), 4 - cluster/ group mem- 
bership, 5 - Stro mgren photometry llStaufferl|priv.com"ml ISone et al.ll200llMSonell2000h . 6 - iNordstromI ||2004|) , 7 - isochrone-fitting (§4.3, 
ISeiss et aklboool) 



Table 5. Sample Statistics 



Category 24 70 fim 24 or 70 fj,m Median Literature 

# Disk # Disk # Disk # F7O/F* Comparison 

excesses Fraction excesses Fraction excesses Fraction total Sensitivity'' 



ALL 


11 


15.7tli% 


8 


ii.4j:^-9% 


15 




70 


2.48 


15%'' 


T. < 5000 K 


5 


18.5+^-^% 


1 


3.7l{-^% 


6 


22.21^-*% 


27 


6.32 




5000 K < T, < 6000 K 


4 




6 


2i.4j:^-5% 


7 


25.0t^'^% 


28 


2.20 


8.5-19%,10-20%'^ 


T. > 6000 K 


2 


13.3+fg3% 


1 


6.7^2^2^% 


2 


13.3tfg3% 


15 


1.74 


33%<i 


single star'^ 


7 


i5.9+ll% 


6 


13.6j:|-»% 


10 


22.7l|-g% 


44 


2.64 


15%=^ 


all binaries 


4 


15.4^8-6% 


2 


7.7+^-^% 


5 


19.21^:8% 


26 


2.29 


23%* 



^The median 24 fim S/N is 57 when including calibration uncertainties. The minimum excess reported is F24/F,=1.09, and consequently the disk 
fractions listed are for F24/F,> 1.1. At 70 ^m, the median sensitivity is determined from the 3-cr excess upper limit for detections and non-detections. 

iLagrange et all ||2000|') 

'^Fractions are at 24 and 70 fim respectively from lMever et aL 1 1 I2OO8I . I2OO5I) 
' ISu et al.l l l2006h 

^Binarity as defined in Table 1; single star — no known companion; spectroscopic or unresolved binary, footnote (f) in Table 1, projected separations 
of 4—78 AU; wide or visual binary - footnotes (e),(g) in Table 1, projected separations of 137-4320 AU. 

'The >3 AU binary separation 70 fim excess fraction (11/48) inferred from Table 4 in jTrilling et alJbOO?!) . The one 24 /im excess source identified 
in lTrilling et al.l ||2007| ) at these separations is also in excess at 70 fim. 



- 60 - 





Fig. 1. — Distribution of detected 24 and 70 fim photometry relative to the expected pho- 
tospheric values for all objects reported in this paper, xro is skewed to positive values since 
we are not sensitive to detect the photospheres for our entire sample, and non-detections 
are not included in the distribution. As a result detections are more likely to be in excess 
than not. The rightmost bin corresponds to the total number of sources with x>10- The 
red vertical dashed line represents the threshold we set for the detection of excess. 



- 61 - 



100 



tin 



tin 



10 



HD 929451 



V 



AU Mic 



1:HD7590 
mHD 73350 



V 
V 



HD 135599;^ 



HD 10008 It 
V HD 59967 

V V ■ 
V HD 1 16956 



HD 17240 



w 



■P HD 1 12429 



V ^ ^ 



V Tf7 



[1 AT Mic 



V V 
V V 
V 

V 



V 



liHD 1835 

V 



HD 186219 

V 



,GJ3400A HD 92139 



Ind 



HD 202730 ™ 1^089 
V HD 177724 



4000 



6000 



8000 



10000 



T, (K) 



Fig. 2. — The observed 70 /iin MIPS flux density (not color-corrected) divided by the ex- 
pected photospheric flux density as a function of effective stellar temperature for our sample. 
3-0" upper limits from non-detections at 70 fim are plotted with blue triangles, and detections 
are plotted with black squares and name labels. The MIPS observations for some sources 
have been previously published (Table 2, Column 8). However, all photometry plotted are 
derived in this work and represent improved absolute flux calibration uncertainties and data 
reduction methods. We are less sensitive to 70 fim excess for cooler stars, introducing a bias 
in our sample. 
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Fig. 3. — The observed K5-[24] color as a function of stellar effective temperature for our 
sample. All colors plotted are re-derived in this work as in Figure 2. The large (>0.2mag) 
uncertainties for several sources are either due to the saturation of the 2MASS images, 
or blending from companions at 24 fim. In green we plot the e mpirical mode l of th e intrinsic 
stellar K5-[24] colors for <5pc main sequence dwarfs derived in lGautier et al.l (120081 ). Sources 
with X24>5 are plotted in blue, indicating possible warm excesses. 



- 63 - 





Fig. 4. — Spectral energy distributions (SEDs) for stars in our sample with 70 fim MIPS 
detections. Stars are ordered from top-left to bottom- right as they are in Table 2. First, 4 
stars with 24 and 70 fim excesses are shown, followed by 4 stars with 70 fim excesses only, 
followed by two stars with 24 fim excesses only and 70 fim detections, and finally stars with 
photospheric 24 and 70 fim detections. Optical and near-infrared photometry are shown 
in black. MIPS color-corrected photometry are shown in blue, and IRAS color-corrected 
photometry are shown in black with 1-a uncertainties shown for detections. Open triangles 
are shown for 3-cr upper-limits. Horizontal error bars on photometry represent effective 
bandwidths. In light brown are the temperature-interpolated Phoenix NextGen sythetic 
spectra fit to the optical and near-IR photometry. In red is the dust component SED, and 
in black is the summation of the synthetic stellar SED and dust component SED depicting 
the total model flux at wavelengths longward of lO/xm. 



- 65 - 




HD 1835 



U K V'\- J H K 



V 



24 tim 70 |J.m 




HD 15089 



U B V i\ ] H K 



24 fim 70 \im 




HD 17240 

U B v'\ J H K 



24 ^ini 70 nm 




HD 92139 



I) B V^c 'i J H K 



24 (im 70 |im 




=1. ^ 



HD 116956 



U B V^c 'c J H K 



V. 



24 urn 70 |im 




(LLm) 



-67- 



1000 - 



100 



< 



10 




3000 



4000 



5000 



6000 



7000 



8000 



9000 10000 



Teff(K) 



Fig. 5. — We plot the distribution of our sources as a function of stellar effective temperature 
and estimated age. Open black circles are stars without MIPS excess; open green squares 
with both 24 and 70 /xm MIPS excesses; open red circles are stars with 70 /xm MIPS ex- 
cesses; open blue diamonds are stars with 24 fim excesses. Uncertainties in age and effective 
temperature are shown as lines centered on each symbol. See §4.3 for discussion. 



Fig. 6. — Left: 24 yum image of TWA 13, with the A and B components labeled in black 
and circled in white. A third background galaxy is also circled. The gray lines divide the 
image into quadrants, and the image is centered on TWA 13A. Right: 70 fim image at the 
expected position of TWA 13. The 70 fim emis sion is coincident wi t h the backg r ound galaxy 
and not TWA 13A as previously published in I Chen et al.l (l2005al ): iLow et al.l (120051). The 
same gray lines are used as in the 24 /im image, centered on the expected position of TWA 
13A. Both the left and right images are to scale, with north up and east to the left (a 20" 
line is shown in black). 



10 100 1000 

Age (Myr) 

Fig. 7. — Fractional infrared excess as a function of estimated age for all stars in our sample. 
Open green squares are stars with both 24 and 70 /xm MIPS excesses; open red circles are 
stars with 70 /xm MIPS excesses; open blue diamonds are stars with 24 /xm excesses; open 
downward triangles are upper limits for stars with no excess derived from the 70 /im flux 
upper limit. Uncertainties in age are shown as horizontal lines. A t;^^ power law decay 
consistent with coUisional disk evolution is shown as a black dotted line. Due to 70 fim 
sensitivities, this plot is incomplete for log(L//j/L=K) < -4.5 and ages less than ~50-100 Myr. 
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Fig. 8. — Fractional infrared excess as a function of measured or estimated stellar rotation 
periods for stars in our sample with T,,<6000 K. We choose 6000 K to be consistent with 
our previous subgrouping of our sample and to select stars with only convective atmospheres 
(there are no stars in our sample between 6000 and 6800 K). Green squares are stars with 24 
and 70 /im MIPS excesses; red circles are stars with 70 fim MIPS excesses; blue diamonds 
are stars with 24 /im excesses. Filled symbols indicate a measured period; open symbols 
indicate a period estimated from v sin i and the radius from the SED fit. 



- 71 - 



s-f 
o 




2 2.5 
Rotation Period (days) / Log (L^/LJ 



3.5 



Fig. 9. — Fractional infrared excess plotted as a function of the negative ratio of the measured 
rotation period to the logarithm of the fractional X-ray luminosity. X-ray activity increases 
to the right, and period increases to the right. Red circles are stars with 70 /xm MIPS 
excesses (HD 135599, HD 7590 and HD 73350); upper hmits for stars without 24 or 70 /xm 
excesses are shown with open black triangles. Only stars with measured rotation periods 
and X-ray fluxes with -5 < Log(Lx/L^,) < -4.3 are included (t^, > ~300 Myr). 



